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Ultraluminous X-ray sources (ULXs) are accreting black holes that may contain the missing population of intermediate mass black 
holes or reflect super-Eddington accretion physics. Ten years of Chandra and XMM-Newton observations of ULXs, integrated by 
multiband studies of their counterparts, have produced a wealth of observational data and phenomenological classifications. We 
review the properties of their host galaxies, list popular spectral models and implications for standard and supercritical accretion 
physics, demonstrate how X-ray timing of these objects places constraints on their masses. We also review multiwavelength studies 
of ULXs, including the optical emission of the binary system and nebulosity around them. We summarize that three classes of black 
holes could power ULXs: normal stellar mass black holes (~ 10 M ), massive stellar black holes (< 100 M Q ), and intermediate 
mass black holes (10 2 — 10 4 M ). We collect evidence for the presence of these three types of compact objects, including caveat of 
each interpretation, and briefly review their formation processes. 
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Q-rl. Introduction 
i 

2 J-l ■ Early identification of a new class of X-ray sources 

An empirical tenet in the study of accreting compact objects 
& is that their maximum radiative luminosity is about Eddington 
— " luminosity, that is the limit where the inward-directed gravita- 
^vq tional force is balanced by the outward-directed radiation force. 
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For spherical accretion of fully ionized hyd rogen, the Edding- 
ton limit can be written as dFrank et alj|2002l) 
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where cr T is the Thomson scattering cross section, m p is the 
proton mass, and M is the black hole (BH) mass. 
| The X-ray luminosity is generally used as a good proxy for 
• * the bolometric luminosity for accreting neutron stars and stel- 
. £h l ar mass BHs. The Eddington argument proved very valuable in 
the early stages of X-ray astrophysic s, to show that most Galac - 
5-H ticQ X-ray binaries are neutron stars dMargon & Ostrikerll973l) . 
Later on, dynamical mass measurements confirmed the exis- 
tence of a population of Galacti c BHs with masses « 5-15M 
dRemillard & McClintockll2006 . for a review): they also satisfy 
the Eddington limit, with peak luminosities < 10 39 erg s _1 . 

The Einstein satellite revealed for the first time a small 
number of non-nuclear X-ray sources in nearby galax- 
ies (particularly, in star-forming galaxies) with appar- 
ent luminosities > 10 39 erg s _1 dLong & van Spevbroeckl 
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19831: IFabbiano & Trinchieril Il987 : Fabbiano Il988l 119891: 
Stocke et alJll991allbh . However, it was not immediately clear 
whether these sources could be a distinct class from Galactic 
stellar mass BH candidates, which were just being discovered 
around the same years. Einstein's relatively crude spatial res- 
olution allowed for the possibility of source confusion in com- 
pact star-forming regions; lack of long-term monitoring made it 
hard to distinguish between persistently luminous sources and 
transient events, e.g., young supernovae (SNe) missed by opti- 
cal observers. 

In the following decade, as nearby galaxies were repeatedly 
observed with X-ray satellites such as ROSAT and ASCA, with 
superior spectral coverage and resolution, it became clear that 
some of those non-nuclear sources could not be SNe, and had 
apparent luminosities well above the Eddington limit of typi- 
cal Galactic stellar mass BHs. For example, the "most likely 
interpretation" of the brightest non-nuclear source in the spiral 
galaxy NGC 1365 was "an ultra-powerful X-ray binary, with 
either a highly super-Eddington low-mass black hole or a very 
massive black hole" with M ~ 100-200M H , "which may pose 
a cha llenge for stellar evolution models" dKomossa & Schulz 
1998h : a third sce nario, beamed em ission, was also proposed 
for other sources (Okad a et al.lll998l) . Thirteen years later, dis- 
tinguishing between those three alternative scenarios remains 
the unsolved fundamental question for this class of sources, 
despite the wealth of new discoveries from Chandra, XMM- 
Newton, Suzaku and Swift, intensive modeling of their X-ray 
spectral and timing properties, and integrating studies of their 
multiband counterparts. 

For a while, different names have been in use to iden- 
tify this class of accreting BHs: extraluminous X-ray binaries 
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Colbert & Mushotzkvl fl 999). superluminous X-ray sources 
Roberts & Warwic 2000), ntermediate-luminosity X-ray ob- 
jects ( Colbert & Ptakl 2002 ). In recent years, the commu- 



nity consensus has settled on ultraluminous X-ray sources 
(ULXs ), adopting the terminology first used by Japanese ASCA 



teams (fokada et al.ll 19981: iMizuno et alll999HMalrishima et al 
2000) 



1.2. Diversity of the ULX population 

A handful ULXs are observed 
X-ray variability at time scales 



to exhibit strong 
as short as minutes 



2007 



( Strohmaver & Mushotzkv 2003; Strohmaver et al 
iHeil et alj|2009l: iFeng et atll2010l: iRao et al.ll2010l) ~nd many 
others have shown rand om variation at long t i me scales from 
days, weeks to years jFeng & Kaaretl 2009; Kaaret & Feng 
2009; iKong et al.1 l20ld; iGrise et al.1 l2010h confirming their 



compact nature. They are 10-100 times more luminous than 
Galactic black hole binaries (BHBs), whose peak luminosity 
is about 10 39 erg s . As the Eddington limit is proportional 
to the BH mass, this suggests that ULXs may contain BHs at 
least 10-100 times more massive than Galactic ones which are 
typically of 10 M , and thus represents a new class of BHs 
with masses in the range of 10 2 — 10 4 M Q , intermediate between 
stellar mass and supermassive BHs. However, some models 
suggest that the Eddington lim it could be viola ted by a factor 
up to ~ 10 for stellar mass BHs (Begelman 2002), and the X-ray 
luminosity may be overestimated due to anisotropic emission 
caused by geome tric beaming expected at high accretion rate 
(Kin g et al.ll2.Q0ll) . In these circumstances, intermediate mass 
black holes (IMBHs) are not required to interpret ULXs. A 
combination of super-Eddington and mildly beamed emission 
from stellar mass BHs could account for ULXs with an 



apparent luminosity up to ~10 41 erg s 1 dPoutanen et alJE007 ). 



Some objects other than BHBs may appear like ULXs if their 
timing or multiwavelength properties are unavailable. Young 
SNe, especially of t ype Iln, could be as lum inous as 10 40 erg s" 1 
in the X-ray band (Im mler & Lewinl l2003). Their X-ray spec- 
trum is usually thermal and soft, and will not show chaotic vari- 
ability. Young X-ray pulsars could theoretically be more lumi- 
nous than 10 39 er g s' 1 but the num ber of such sources appears 



to be very small dPerna et al.l 120081) . In general, these objects 
can be distin guished from accreting sources with variability 
(Kaaret 2008), but it is possible that a SN leaves behind a com- 
pact remnant which becomes X-ray luminous when a fallback 
disk forms (|Li 120031) and mixed properties may be seen. In a 
survey of ULXs in nearby galaxies with archival XMM-Newton 
data, it was found that the X-ray spectra of a few sources can 
be decomposed into a power-law component plus a thermal 
plasma, and in one of them X-ray variability was detected, sug- 
gesting th at they may be recent SNe with an accreting com- 
pact core dFeng & Kaaretl 120051) . Background active galactic 
nuclei (AGN) that happen to lie behind a nearby galaxy could 
be misidentified as ULXs in the galaxy. Their cosmological 
location c ould be revealed if redshift is measured from optical 
lines (e.g. .iFoschini et al.ll2002bt iMasetti et ai1l2003l: lArp et al.l 
2004i lClarketal.ll2005l) . The fraction of backgr ound objects is 



2005). Many ULXs are found to be spatially associated with 
nebulae, young star clusters, or globular clusters; these are most 
likely located in the galaxy rather than foreground or back- 
ground objects. Even if they are identified as BHBs, their nature 
is still not unique in terms of the companion star (low mass or 
high mass), the mode of the mass transfer (Roche-lobe overflow 
or wind-fed), or the mass of the accretor (stellar ma ss or inter- 
mediate mass). For example. IFeng & Kaaret ( 2006a ) found that 
ULXs in young, star forming galaxies are much more variable 
than in old, elliptical galaxies (associated with globular clus- 
ters). 

No matter what their nature is, ULXs are of great interest 
in astrophysics. If they are IMBHs, their formation mechanism 
could be rather different from typical stellar mass BHs; extreme 
environment, more massive and/or less abundant progenitors 
are lik ely required ( Mad au & Reesll2001 ; Portegies Zwart et al. 
2004). IMBHs could have played an important role in the 



usually high in massive, elliptical galaxies (e.g., Tang & Zhang 



formation of supermassive BHs in the early universe and be 
important targets to search for gravi t ational wave emission 
( lEbisuzaki et alj|200ll; I\^lonterill2010t: lMillerll2002b . Instead, 
if they are stellar mass BHs, it is interesting to investigate why 
they are so distinct from their Galactic cousins. They may be 
good candidates for the study of t he accretion physics at near or 
super Eddington accretion rates (lOhsuga et al. 2005), and they 
may shed light on binary stellar evolution channels not seen in 
our Galaxy. 

1.3. Definition of ULXs 

ULXs are usually defined as non-nuclear, point-like objects 
which at least once have been observed at an apparent isotropic 
X-ray luminosity higher than that of stellar mass Galactic BHs: 
typically, Lx > 10 39 erg s _1 in the 0.3-10 keV band. Some au- 
thors (e.g. Kaaretl 12008) prefer to use Lx > 3 x 10 39 ergs s _1 , 
that is the Eddington limit of t he heaviest stellar BH (20M H ) 
expected at normal metallicity ( Belczvnski et al. 2010h . to ex- 
clude a large number of "faint" ULXs that should be similar to 
normal Galactic BHBs. 

Such an empirical definition may have ambiguities and also 
include young X-ray pulsars, SNe, and SN remnants as men- 
tioned above. The consensus in the ULX community is to re- 
strict the definition of ULXs to accreting BHs. Therefore, we 
define ULXs as non-nuclear accreting BHs with peak luminosi- 
ties inferred assuming isotropic emission above the Eddington 
limit of normal stellar mass BHs. This definition is confined 
to compact objects that are most likely accreting from a com- 
panion star in a binary system, but more exotic possibilities, for 
example isolated IMBHs accreting via the Bondi process from 
the surrounding medium are not excluded a priori. Observa- 
tionaly, variability information is required to define a ULX. 

Recently, many ULXs are f ound to have a mechanic al power 
output as high as in radiation dPakull & Mi rioni 2003). In one 
of them, the mechanical output is domina nt while the radia - 
tive output is way below the ULX threshold dPakull et al.l2010h . 
This source is apparently not a ULX, but speculated to be an off- 
set ULX or quiescent ULX (see SectionQfor details). We thus 
define them as "ultrapowerful sources (UPSs)" for complete- 
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ness, to standard for the same population of sources perhaps in 
a different state or viewing angle. 

1.4. ULX Catalogs 

A ROSAT/HRl catalog of luminou s, non-nuclear X-ray 
sources bv lColbert & Mushotzkv ( 1999) proved to be very im- 
portant for bringing ULXs to the forefront of research, and for 
emphasizing the possible observational link with theoretically- 
predicted intermediate mass BHs, with masses 10 2 < M < 
1O 4 M . Larger and more systematic X-ray source cata- 
logs based on re-analysis of archival ROSAT/HR l data soon 
followed: in particul ar, by iRoberts & Warwick] ( |2 000). by 
Colbert & PtakN2002h . and bv lliu & Bregmanl d2005l) ~which 



included 106 sources above 10 39 erg s 



The first XMM- 
A Chan- 



Newton catalog was done by Foschini et al. (2002a). 
dra catalog with 15 4 sources in 82 galaxies was produced by 
Swartz et al. ( 20041). and 226 objects were included in the cata- 
log by Liu & Mirabel d2005h . compiled from mixed literature 
sources. Such works allowed the first population studies of 
these objects as a well recognized group, and in relation to the 
properties of their host galaxies. They also provided the ratio- 
nale for longer, targeted observatio ns with Chandra a nd XMM- 
Newton. Today, the largest catalog (Walt on et al ■1201 ll) is based 
on the 2XMM Serendipitous Survey, and contains 470 ULXs, 
of which 367 not listed in previous compilations. Updated cat- 
alogs based on Ch andr a observati ons have been compiled by 
Swartz et al.1 d201ll) andlLiul d201lb . 



2. Emission Scenarios and Black Hole Masses 

2.1. Different ways of making a ULX 

There are three basic ways to explain non-nuclear sys- 
tems with apparent accretion luminosities > 10 39 erg s _1 , 
using only standard accr etion physics already included in 
Shakura & Sunvaevl(ll973l) , s model. We can either increase the 
mass of the BH (thus, increasing its Eddington limit), or at- 
tribute the enhanced brightness to beamed emission, or require 
a mass accretion rate high enough to allow the source to exceed 
the Eddington limit by a factor of a few. It is also possible that 
the most luminous ULXs may be explained by a combination 
of all three factors. Assuming solar abundances, the apparent 
luminosit}0 of an accreting BH is: 
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dShakura & Sunvaevl[l973l: IPoutanen et al]l2007l) . where b < 1 
is the beaming factor, and m is the dimensionless accretion rate 
at large radii (typically, the mass transfer rate from the donor 
star), that is normalized to the Eddington accretion rat^§ m 



\ bat is, the observed flux x4nx distance . 
3 This definition implicitly includes a "standard" radiative efficiency r\ = 0. 1, 
typical of efficient disk accretion; hence, a dimensionless accretion rate m ~ 1 
produces a luminosity L ~ Z-Edd- Other authors prefer to define m = Mc 2 /Leh, 
in which case the Eddington limit is reached at m ~ 10. 



M/MEdd = O.lMc 2 /LEdd- Let us examine the main tenets of 
these alternative or complementary scenarios. 

Strong beaming (1/b » 1) — A decade ago, ac- 
creting BHs with strongly beamed emission were already 
well kn own in extragalactic astronomy (BL Lac objects and 
blazars: lUrrv & Shaferlll984l) . and at least one Galactic BH 
(SS 433) was known to have semi-relativistic jets, massive 
outflows and likely anisotropic emission. Naturally, a pos- 
sible explanation for the early sample of ROSAT ULXs was 
that they were microblazars, or SS 433-like objects seen 



face- on (Fabrika & Mes chervakovll2001t iKording et alJ 12002: 



iBegelman et alJl2006l) . The relativistic beaming scenario soon 
ran into several fatal difficulties. It predicts that for every ULXs, 
there should also be a large number of lower-luminosity beamed 
sources: « 30 sources with apparent L ~ 10 39 erg s _1 for ev- 
ery ULX at L ~ 10 40 erg s" 1 ; but ther e is no evidence of this 
population ( Davis & Mushotzkv 2004). Instead, empirical lu- 
minosity functions show only * 5-10 sources m ore luminous 
than 10 3 9 erg s" 1 for eve r y ULX at 10 40 erg s" 1 dWalton et al 



2011; Sw artz et al.1 2004; Grimm et al. 2003), consistent with 



the higher end of the high-mass X-ray binary luminosity func- 
tion. There is also no multiband evidence of the large number 
of jet sources beamed in other directions, required by this sce- 
nario (Dav is et al .120031) . Conversely, the presence of photoion- 
ized bubbles around several ULXs (for example Holmberg II 



X-l: 



Pakull & Mirionil2002l;lKaaret et alj2004l:lLehmann et al 



l2Q05t NGC5408 X-l: Kaaret & Corbelll2009l) requires quasi 
isotropic X-ray emission with luminosities « 10 40 erg s _1 . Fi- 
nally, the majority of ULXs do not have radio counterparts 
and/or rapid X-ray variability, expected for relativistic beam- 
ing. For all those reasons, relativistic beaming is now ruled out 
as a general scenario. 

Mild beaming and/or super-Eddington accretion (l/b < 10, 
m » 1) — When the accretion rates exceeds the Edding- 
ton limit, radiatively-driven outflows are launched from the in- 
ner part of the disk, inside the sp herization radius where the 
disk becomes geometrically thic k d Shakura & Sunyaev 1973 
| Poutanen etly 20071 iKiiT'' 



y thic k Qbhakura & sunvaev IV /i 
2009). It was suggested (Ki ng et al 



20011: lKindl2009h that the wind may produce a funnel wall 



that scatters and collimates the emission along the axis per- 
pendicular to the disk plane. The apparent high luminosity 
of ULXs is th en a comb ination of this mild collimation (with 
b ~ 70/m 2 . lKind r2009) and mild super-Eddington emission 
(L ~ LEdd(l + 3/5 lnm)). In this scenario, accretion rates 
m ~ 10-30 are enough to explain most or all ULXs up to lumi- 
nosities as 10 41 erg s~' with BH masses < 2OM . 

Recent n umerical r ad iation -m agneto-hydrody n amica l 



simulations (|Ohsuga et alJ 120091: Takeuchi et al. 2010; 
Mine shige & Ohsugal 1201 11) supporte d and advanced earlier 
analytical results (slim disk models; lA bramowicz et al. 1988; 
Wata rai etal1l200U lEbisawa et all 120031: iHeinzeller & Duschl 
120071) . and t he radiation-hydrod ynamical simulations of, 
among others, lOhsuga et alJ d2005l) . An important result of 
the simulations is that even a moderately super-critical mass 
supply m w 5 can produce a total luminosity m 1.7LEdd (in 
rough agreement with the standard analytical scaling used 
earlier), and an apparent luminosity w 22Lsdd for face-on 
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observers, due to mild collimation in the magnetized o utflow 
dMineshige & Ohsugall201 lHOhsuga & Mineshigell201 ll) . 

An alternative explanation for ULXs relies on super-critical 
accretion but does not require collimation of the emerging ra- 
diation. It was suggested that radiation pressure-dominated ac- 
cretion disks have strong density inhomoge neities due to th e 
development of photon-bubble instabilities dBegelmanl 12002 ) . 
Radiation would escape from such disk regions at rates ~ 10 
times above the classical Eddington limit for stellar mass BHs. 

Quasi-isotropic Eddington luminosity (1 jb > L^Edd ~ 1) — 
Although some models and simulations offer us compelling rea- 
sons to ignore the Eddington limit as a constraint for the ULX 
interpretation, some authors remain wary of this scenario. They 
point out that the Eddington limit works well for Galactic BHs 
and neutron stars. More significantly, the isotropic bolometric 
luminosity of the whole population of > 60, 000 quasars at red- 
shift 0.2 < z < 4 in the Sloan Digitized Sky Survey is strongly 
Eddington-limited at all BH masses, with only a very small 



masses » 5-15M (for reviews see McClintock & R emillard 



fraction of sources between 1 and 3 LecM (Steinhardt & Elvis 



20101) . Perhaps like the bumblebee, the Eddington limit should 
not fly, but does. 

Thus, it is argued that the null hypothesis for ULXs should 
be that their luminosities are » LEdd, despite the difference on 
the gas reservoir between AGN and binary systems, and their 
radiative anisotropy is only a factor of > 1 (a standard accre- 
tion disk has a beaming factor l/b - 2 in the absence of any 
other collimation). These conditions are possible to explain 
most ULXs up to a few xlO 40 erg s _1 with accreting BHs of 
masses < 100M o , which is the theoretical upper limit for a 
direct collapse of a single stellar core in the current universe 
(see Section 12.2b . Sources more luminous than that may re- 
quire IMBHs. For example, HLX-1 in ES0243-49, the most 
luminous ULX found to date, has a peak bolometric l uminosity 



10 42 erg s" 1 (IFarrell et al J 120091: Davis etalJl201lh . If it is 



Eddington limited, a black hole mass of at least 1O 4 M is in- 
ferred. Even if super-Eddington emission is allowed, an IMBH 
is required. 

Quasi-isotropic sub-Eddington luminosity (l/b ~ 1, m < 
1) — Finally, it is possible that ULXs are accreting at sub- 
Eddington accretion rates with nearly isotropic emission or disk 
beaming, which requires the compact objects to be IMBHs 
(1O 2 M < M < 1O 4 M ). Then we expect to see similar state 
transition and X-ray spectral properties from both ULXs and 
Galactic BHBs, scaled only by the mass. 

2.2. Three types of non-nuclear black holes 

The emission scenarios outlined above require three different 
ranges of BH masses, with different formation mechanisms. 

Ordinary stellar mass BHs (M < 2OM ) — The first ad- 
vantage of explaining ULXs with ordinary stellar mass black 
holes (sMBHsS) is that we already have direct evidence of 
their existence in the local universe. There are more than 20 
BHs identified in X-ray binaries with dynamically measured 



2006^ iRemillard & McClintockll2006l) . In the Milky Way, the 
most massive sMBHs are those in GR S 1915+105 (M 
1 4.0 + 4.4 M ; lHarlaftis & Greinerll2004 and in Cygnus X-l 
(M = 14.8+1.0 M^: IOrosz et al.l201 ll) . A si milarly massive BH 
is found in M33 X-7 (M = 15.65 + 1 .45 M dOrosz et~alll2007t) . 
An a ccreting binary system in the nearby starburst gala xy IC 
10 dPrestwich et alj|2007t ISilverman & Filippenko]|2008l) con- 
tains a BH with a mass of « 2O-4OM (the mass uncertainty is 
mainly due to the unconstrained companion mass and inclina- 
tion angle of the binary plane). The second advantage is that 
they are the easiest BHs to form: standard stellar models pre- 
dict sMBHs as the endpoint of individual massive star evolu- 
tion, for a wide range of stellar masses and metallicities. Thus, 
the strictest constraints to sMBH models for ULXs come not 
from the BH formation process, but from the extreme super- 
Eddington accretion, beaming and outflows required by these 
models. 

Massive stellar BHs (2OM < M < 1OOM ) — The mass 
of a compact remnant depends both on the initial mass of the 
stellar progenitor, and on how much mass is lost through stel - 
lar winds during the stellar lifetime (iFrver & Kalogerall200ll) . 
Radiatively-driven winds from massive stars are strongly de - 



pendent on their metallicity (IVinkl 1201 It IVink et al.l 1201 ll) . 



Therefore, metallicity plays an important role in stellar evo- 
lution and determines the maximal mass o f a BH that can be 
formed via core collapse. Belczvnski et al. ( 2010h argued that 
the collapse of a single star can form a BH up to as 15M at 
solar metallicity, « 3OM at 0.3 times solar metallicity, and 
~ 8OM at 0.01 times solar metallicity; consistent results (upper 
BH mass limit * 7OM at low metallicities ) that a re obtained 



4 To avoid confusion with the acronym "SMBH" usually adopted for super- 
massive black holes. 



bv iHeger et al.l (120031) and lYungelson et al.l (120081) . These re- 
sults can well explain the inferred BH masses in Galactic X-ray 
binaries (average metallicity close to the solar value) and in the 
IC 10 source (metallicity as 0.3 solar). In summary, current stel- 
lar evolution models predict a BH remnant distribution peaked 
around 1OM , with a tail extending up to < 1OOM . As the 
formation of stellar BHs more massive than w 2OM may re- 
quire rarer, more extreme conditions, we will refer to them as 
"massive stellar black holes (MsBHs)". 

There is another interesting physical difference between 
sMBH and MsBH formation, in addition to different progenitor 
mass and metallicity: sMBHs are probably formed through a 
SN explosion and subsequent rapid fallback; MsBHs are prob- 
ably formed from direct core collapse without an explosion 
(!Frverjll999l) . This is because for the more massive objects, 
the gravitational energy released during core collapse is not 
enough to unbind and expel the stellar envelope. At low metal- 
licities, th e threshold for direct collapse is an initial stellar mass 
~ 40M Q dFrverlll999l: IHeger et alJl2003h : at solar metallicity, 
even the most massive stars may lose enough mass in winds 
and explode as SNe. 

Intermediate mass BHs (M ~ 10 2 -10 4 M ) — To produce an 
IMBH of ~ 1O 2 -1O 4 M , the core collapse of an isolated star in 
the current epoch is not a viable process. This is not because of 
the lack of very massive stars (in fact, several s tars with initial 
masses ~ 150-300M o may exist in the LMC; Crowther et al 
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2010), but because Helium cores more massive than a 70M o 
(corresponding to initial stellar masses > 14OM at zero metal- 
licities) undergo an electron/positron pair-instability explosion 
that lea ds to the complete disruption of the star, w ithout BH for- 
mation dBond et alJll984tlHeger & Wooslevll2002l) . However, a 
BH may form again if the Helium core is more massive than 
« 13OM , corre sponding to initial stella r masses > 26OM at 
zero metallicity dHeger & WooslevlEoolh . 

It has been suggested that metal-free population III stars 
formed in th e very early U niverse could reach masses of a few 
hundred M Q ( Larsonll998 ). above the pair-instability limit, and 
thus may have collapsed into IMBHs dMadau & ReeslEoOl ). In 
young and dense star clusters, dynamical friction could lead to 
massive stars sinking towards the center and undergoing run- 
away collisions and mergers on timescales < 10 6 yr. Numer- 
ical simulations suggest that such a process may result in the 
formation o f a huge star, which may subsequently collapse into 
an IMBH dPortegies Zwart & McMillan! l2002t iGurkan et al.1 
l2004tlPortegies Zwart et al.l2004IVanbeveren et al.l2009t) . Dv^ 
namical evidence for IMBHs with masses ~ 10 4 M o has 
been proposed for a few globular clusters, e.g., Gl in M31 
(Gebh ardtetaill 2005). although t h e issue is still controv ersial 
dAnderson & van der MarefeoiOh . IVesperini etail d2010h sug- 
gested that in massive and/or compact globular clusters, a cen- 
tral seed sMBH may grow by up to a factor of 100 via accre- 
tion of gas lost by the first generation of cluster stars in their 
red-giant phase. IMBHs may also wander in the halo of major 
galaxies, after tidal stri pping of merging sate l lite dwarfs that 
conta ined nuclear BHs (IKing & Dehnenll2005 ; Bellovarv et al 



2010h : this is currently the most plausible explanations for the 



most extreme ULX found to date, HLX-1 in ES0243-49. 

To summarize, three types of non-nuclear BHs could be pow- 
ering ULXs (Table [1). 

Table 1 : Three types of non-nuclear BHs that could power ULXs by accretion 



Type 


Mass 


Progenitor 


l/b-L/L^ 


sMBH 


-10 M 


normal star 


> 1 


MsBH 


<100M o 


low-Z star 


~ 1 


IMBH 


10 2 -10 4 M a 


pop III 
cluster core 
stripped nucleus 


< 1 



3. Population studies 

3.1. Host galaxy environment 

Assuming that ULXs are binary systems with a BH receiv- 
ing matter directly from an individual donor star, we can try 
to fit them into the traditional classification of Galactic X- 
ray sources, as either high-mass or low-mass X-ray binaries 
(HMXBs and LMXBs, respectively). HMXBs have a massive 
OB donor star with a characteristic age < 10 7 yrs; LMXBs 
have a low-mass donor (often in the subgiant or giant phase) 



with ages > 10 9 yrs (e.g. iFrank et al.1 l2002h . HMXBs are 
mostly found in spiral and Irr galaxies, with current or re- 
cent star formation, while LMXBs populate old spheroidals and 
globular clusters. ULXs are found both in ellipticals and in 
spiral/irregular galaxies; however, the most luminous ones are 
found predominantly in star-forming galaxies. 

Specifically, more than two thirds of ULXs found in ellip- 
ticals ar e in the lowe st luminosity range, with Lx < 2 x 10 39 
erg s _1 (ISwartz et al.ll2.Q04l) , and can easily be explained as the 



upper end of the LMXB distribution, with BH m asses ~ 2OM . 
It was even suggested (llrwin et al. 2003, 12004 that all ULX 



candidates in ellipticals with apparent Lx > 2 x 10 erg s 
could be dismissed as background or foreground sources. To- 
day, it is generally accepted that there are a few genuine ULXs 
in ellipticals at luminosities ~ 2-10 x 10 39 erg s _1 , but none 
has been found above 10 40 erg s -1 . In spirals, by comparison, 
one third of ULXs have luminosities > 4-5 x 10 39 erg s~\ and 
about 10% ha ve luminosities > 10 40 erg s _1 dSwartz et al . 2004; 



Walt on etalj|2011l) . Both early-type (Hubble types to 4) and 
late-type (Hubble types 5 to 9) spirals may contain ULXs, with 
similar proba bility, w i thin t he current statistical error of the 
samples: I Walton et alJ ( 201 1 ) fi nd slig htly more ULXs in early- 



type spirals, while Swar tz et al.l (1201 II) find slightly more in late 
types, due to the different spatial resolutions of the instruments. 

Furthermore, the fact that a luminous ULX is found in a spi- 
ral galaxy does not imply that it has an OB donor star. Many 
ULXs in spirals are not located in regions (within a few 100s 
of pc) with eviden ce of current or very recent star formation 
(Mushotzky 2006). Using the optical color distribution in a 
sample of face-on host galaxies. ISwartz et al. (2009) found that 
the brightest ULXs turn on ~ 10-20 Myr after the end of star 
formation. Characteristic ages > 10 Myr for the surrounding 
stellar population were al so found for example in the case o f 
NGC 1313 X-l and X-2 (lYang etalJIioill iGrise et alJl2008h. 
IC 342 X-l dFeng & Kaaretll2008l). NG C 4559 X-l dSoria et al 
20051). and NGC 1073 IXO 5 dKaaredEooll) . 



Swartz et al 



(2011) normalized the radial distance of all 
ULX candidates in their sample to the deprojected D25 diam- 
eter of their host galaxies. They found that the ULX surface 
density profile is strongly peaked towards the center. The pro- 
file can be fitted with a Sersic function, with a best-fit index 
consistent with a De Vaucouleurs profile (n = 4), and an expo- 
nential scale-radius ~ 10~ 3 Z)25. This finding suggests that the 
ULX population is intrinsically related to the stellar population 
of their host galaxies. Less centrally concentrated distributions 
would be expected if most ULXs were primordial halo relics 
such as Population III remnants, or nuclear BHs of accreted 
satellite galaxies. 

A Chandra sample of ULXs in 58 face-on star-forming 
galaxi es within 15 Mpc, studied in details by Swartz et al.l 
d2009l) . does not show statistically significant associations be- 
tween ULXs and young, massive (M» > 1O 5 M ) star clus- 
ters, with in 100" x 100" regions. The same conclusion was 



found by Zezas & Fabbiano (2002) for ULXs in the Antennae 
galaxies. This is inconsistent with formation scenarios based 
on runaway core collapse and mergers of O stars inside a com- 
pact, young stellar cluster, leading to the possible formation of 
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of the resulting superstar ( 


Porteeies Zwart & McMillan 2002; 


Porteeies Zwart et al. 2004 


; Freitag et al.1 2006|). In fact, most 



ULX host galaxies do not even have stellar clusters sufficiently 
massive and compact to satisfy the requirements for runaway 
core collapse. In theory, IMBHs could be formed inside clusters 
that have since dispersed. However, the evaporation timescale 
of such clusters would be too long to explain the observed as- 
sociation of many ULXs with young (< 20 Myr) stellar popu- 
lations. 



3.2. Specific frequency 

Here, we further quantify the asso ciation between ULX pop- 
ulations and host galaxy properties. ISwartz et al.l (120081 1201 II) 
studied 107 ULXs in a complete sample of 127 nearby galaxies 
(d < 14.5 Mpc), selected to be above the completeness limit of 
both the Uppsala Galaxy Catalog and the Infrared Astronomi- 
cal Satellite (IRAS) survey. The total star formation rate (SFR) 
of the galaxies in the sample is ss 5QM G yr _1 based on the 
IRAS data, consistent with other measure ments of star forma - 
tion based on the Ha and [O n] indicators dSwartzetalJl201lh . 
This, in turn, corresponds to an SFR density « 0.01M G yr~' 
Mpc 3 in the local universe. Finding a 110 ULXs in that vol- 
ume, after correcting for incompleteness and background inter- 
lopers, implies that we expect, on average, » 2 ULXs for a SFR 
of IM Q yr" 1 , and 1 ULX in a volume of * 50 Mpc 3 . This 
is consistent with the HMXB lu minosi ty function in local star- 
forming galaxies: iGrimm et al.l Q2003) find w 30 sources with 

39 



L2-10 ^ 5.5 x 10 erg s (correspo nding to Lqj iq > 10 



erg s in the ULX definition used in ISwartz et al.M2.01 II) in a 
sample of galaxies with a total SFR ss 16M yr -1 . 

The average specific frequency in star-forming galaxies, in 
the local universe, is 1 ULX per a* 1O 1O M , for galaxies 



with dynamical m ass of 1O 1U M (Sw artzet al.l I2008L 12011 



Walton et al.ll2011l) . T his rate decreases for increasing galax y 
mass, scaling as M" 06 dSwartz et alj|200c^lwaiton et alj|201 lb . 



In elliptical galaxie s, the specific rate is 1 ULX per * 10 n M o 
dWalton et alHoI 1|). Thus, it appears that low-mass spirals and 
irregulars are more efficient at producing ULXs, per unit mass. 
This is true at least down to galaxy masses ~ 1O 85 M . 

There may be two reason s why low-mass ga laxies contain 
more ULXs per unit mass dWalton et al.l 1201 II) . One is that 
their SFR per unit mass is higher, so they are more efficient 
at producing ULXs. The other is that they have lower metal 
abundances, which may favor the formati on of heavier stellar 
BHs from the collapse of massive O stars dZ ampieri & Roberts 



2009; iMapelli etalj|2009t ISoriaet all 12005 



Pakull & Mirioni 



2002). However, it is suggested by Prestwich et al. ( 2010h that 



the metallicity does not affect the ULX production strongly af- 
ter comparing two pairs of colliding galaxies. 

Finally, it is often noted that many luminous ULXs reside 
in colliding or strongly interacting galaxie s: the Antennae 
dZezas et al] 120021: IZezas & Fabbian d l2002h. the Cartwheel 



dGao et al.ll2003l:IWolter & TrincMeriKOCMt IWolter et alJbOOd 
ICrivellari et alJ l2009h . the Mice dReadl l2003h. NGC 4485/90 



dRoberts et al 



2002), NGC 7714/15 dSoriaetalJ 12004 



— ^^^^^ 








^ M = 103.9 L 39 


0.8 




N = 97.6 L 3a -<>- 5 exp(-L 39 /1 9.7) \L 











SmithglalJ l2Q05|), NGC 3256 dLira et alJ l2002b . and a few 



1 2 5 10 20 

L(0.3-10keV)/(10 39 ergs- 1 ) 

Fi gure 1: Cumulative ULX luminosity function from the Chandra survey 
of Swartz e t all feOllh . For this sample, a cut-off power-law with a break 
a 2 X 10 40 erg s _1 matches the observed distribution better than an unb roken 
power-law, although the issue remains controversial (cf. Walton et al. 201l|). At 
the low-energy end, it also mat ches the slope and normalization of the HMXB 
distribution below 10 39 erg s _1 jGrimm et alj2003l) . 



more Arp galaxies. This may simply be due to the fact that 
colliding galaxies have the largest SFR, and the number of 
ULXs scales with the recent SFR. 

3.3. Luminosity functions 

The cumulative luminosity function (LF) of a population of 
accreting sources (also known as logjV(> S)-logS curve) is 
expected to have a break or cut-off at the Eddington luminosity 
of the most massive objects in the population. Therefore, in the 
absence of direct mass measurements, the LF can be used as a 
rough indicator of the maximum mass. In fact, things are not 
so simple. The LF of X-ray binaries in our Galaxy and in other 
nearby galaxies does not show any features at w 2 x 10 38 erg 
s _1 , Eddington limit of the neutron star population. Moreover, 
the location of the upper cutoff may be a function both of mass 
(via the Eddington limit) and of the age of the population, s ince 
the most luminous sourc es tend to have a shorter lifetime dWu 
200 It iKilgard et alj|2002l) . Thus, careful modeling is required 
to interpret a LF. 

With this caveat, we can use the ULX LF to test the sce- 
nario that ULXs in star-forming galaxies are mostly the high- 
luminosity tail of HMXBs (plus some additional LMXBs), and 
ULXs in ellipticals are the high-luminosity tail of LMXBs. 
Several Chandra and XMM-Newton surveys of discrete X- 
ray sources in nearby galaxies have prod uced consistent re- 
sults in agreement with this interpretat ion (Grimm et al. 
Colbert et : al]l2004l: ISwartz et alJl2004l: IWalton et alJl201 ll) . In 
summary, they have all found that the cumulative LF of young 
X-ray populations in star-forming galaxies is a simple power- 
law with a slope * 0.6-0.8, while the cumulative LF of old 
X-ray populations in ellipticals and spiral bulges has a slope 
of ~ 1.5. Thus, short-lived HMXB-like sources dominate the 
high-luminosity end of the distribution. The normalization of 
the HMXB LF is proportional to the SFR; the normalization of 
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tne lmab J_,r is propc 
2004lGirfanovll2004 . 



the LMXB LF is prop ortional to the stellar mass (IGilfanov et al.l 4.1. Power-law spectra: hard or steep power-law state 



What remains controversial is the behavior of the LF for the 
young ULX population above a n X-ray luminosity » 10 40 erg 
s _! . From their Chandra survey. ISwartz et al.l (1201 II) find that a 



pure power-law N(> Lx) ~ L x a provides a very poor fit to the 
data above that threshold (Figure^- A cut-off power-law model 
with N(> Lx) ~ Ly? exp(-Lx/L c ) provides an acceptable fit, 
with a « 0.5 (as required by the low-luminosity end of the LF) 
an d L c » 2 x 10 40 erg s " 1 . A very sim i lar re sult was also found 



by iGrimm et al.l d2003l) . ISwartz et all d2004 . and iMineo et al 
The cut-off power-law model does not preclude the 



possibility of finding ULXs above 2 x 10 erg s , but th ey 



would be increasingly rarer. It predicts (ISwartz et al.ll201 ll) 1 
ULX with L x > 5 x 10 4() erg s 1 within a radius of ~ 13- 
29 Mpc, and 1 ULX with L x > 10 41 erg s -1 within a radius 
of ~ 42-119 Mpc. If the lumino sity of a few ULXs recently 
discovered by ISutton et al.1 (1201 ll) is confirmed to be w 10 41 
erg s 



, their existence would not be consistent with the cut- 



off model of ISwartz et al .1 (1201 ll) . and may suggest a p hysically 
different population of accretors dKing & D ehnen 2005). And, 
of course, so does HLX-1 in ES Q243-49 with i t s peak out- 
burst luminosity « 10 42 erg s _1 dFarrell et al . 2009; Davis et al 
201 lh. 



In contrast with t hose findings, the XMM-Newton survey of 
Walt on et al.l (1201 ll) does not see any turnover or cut-off in the 



cumulative LF above 10 40 erg s , but their data cannot rule 
out the presence of a similar break. At this stage, the dis- 
crepant result may be caused by a diff erent binning techn ique, 
or by small-number statistics. Both fhe lWalton etal. ( 201 1 ) and 
Swartz et al.l d201 lb complete subsamples have only 5 sources 



with L x > 2 x 10 40 erg s" 1 out of a total of * 100 ULXs. A few 
other ULXs above that threshold are known, but cannot simply 
be added to the LF because they do not form part of a complete 
sample. 



4. X-ray Spectral Properties 

Thanks to the large effective area of XMM-Newton, combined 
with Chandra observations, the spectral and timing properties 
of ULXs in the 0.3-10 keV band have been widely investigated 
and modeled over the last decade, even though the physical in- 
terpretation is still controversial. Most of the X-ray studies have 
focused on bright ULXs (fx > 10~ 12 erg crrT 2 s _1 ) in nearby 
galaxies (< 10 Mpc) that allows spectral and timing analysis 
at high signal to noise ratios. In the next two sections, we re- 
view X-ray spectral and timing properties of bright ULXs and 
discuss their physical interpretation in comparisons with well- 
defined emission properties of Galactic BHBs. 

Phenomenologically, the energy spectra of ULXs in the 0.3- 
10 keV band can be divided into two groups: those that are 
consisten t with a simple p ower-law, and those that are more 
complex (Ma kishimair2007l) . The latter mainly refers to a mild 
broad curvature (convex shape) over the whole band, a break or 
steepening above ~ 2 keV, or a soft excess below ~ 2 keV. A 
typical ULX spectral shape of the latter is sketched in Figure[2] 



Many ULX spectra are well fitted by a single, absorbed 
power-law model. The pho ton index has a broad distribu- 
tion, peaked at T » 1.8-2 .0 dSwartz et all 12004 Iwinter et al 



2006; iBerghea et al. 2008) but with a few much harder (down 
to T « 1) as well as much softer sources (up to T * 3). 
The distribution does not appear to be bimodal; there is no 
evidence of a gap b etween a "hard state" and a "soft state". 
Berg hea et al.l d2008l) found that the power-law photon index 



tends to be harder at increasing luminosity; a selection of ULXs 
with high luminosity (Lx > 10 40 erg s ' 1 ) an d hard power-law 
spectrum (T < 1.8) is also presented in lSorial (1201 ll) . 

Some of the hard sources with multiple observations show 
strong flux variability (so metimes by an order o f magnitude) but 
constant photon index dFeng & Kaaretl 120091; [Kaaret & Feng 
l2009t ISoria & Ghoshl |2009|). Similar beh avior has bee n seen 
from Galactic BHBs in the h ard state dBelloni et al.l 12005 : 
iRemillard & McClintockll2006l) . The X-ray lumi nosity in the 
hard state is usually < 0.03 times Eddington dFender et al 



2004), but during hard-to-soft state transitions t he hard state can 



reach a luminosity up to 30% of Eddington (Miva kawa et al 



|2008HYu & Yaiill2009l) . Thus, if hard power-law ULXs belong 
to the same canonical hard sta te defined in Galac tic BHs, their 
BH masses must be > 1O 3 M d Winter et al.ll2006l) . 

A key signature of the hard state is the presence of com- 
pact, continuous radio jets. So far, the only possible detec- 
tion of such jets is from IC 342 X-l, with VLA observations 
( Cseh et all 20 1 1 ) . Usi ng the fundamental plane of accreting 



Cseh et aiT ( 201 1 ) estimate a compact object mass (1.2- 



BHs, 

13.6) x 1O 3 M . Further observations are needed to spectrally 
constrain the nature of the radio core in this system. 

M82 X-l is another good candidate for hard state ULXs. 
This sour ce displays low frequency qu asi-periodic oscillations 
(QPOs) dStrohmaver & M ushotzkv 20Q2) above a band-limited 



broad noise component, and its X-ray sp ectrum is a feature - 
less power-law with a photon index « 1 .7 dKaaret et al.ll20 06a). 
We note that M82 X-l may be one of the very few ULXs that 
change their spectral state during outbursts, switching from a 
hard to a thermal state (see Section l4~3l >. 

The hard spec trum may come f rom an advection-dominated 



accretion flow (E sin et al.l 11997). a luminous hot accretion 



flow dYuanl |20 - 01 |). a classica l hot corona above the accre- 
tion disk dLiang & Price 1977 ), or from a combination of syn- 
chrotron and invers e-Compton emission at the base of the jet 
dMarkoff et al.ll200"l ). The first mechanism is limited to low ac- 
cretion rates, but the latters can in principle opera te also near the 
Eddi ngton luminosity and be relevant to ULXs dFreeland et al 



2006). Thus, an alternative interpretation favored by some au- 
thors is that hard power-law ULXs are not in the canonical 
low/hard state; in other words, that the ULX hard state is not 
confined to luminosities < 0.1 Eddington, and does not switch 
to a softe r, disk-domi nated thermal state as the accretion rate 
increases (Sorial l201 ll) . We could call it a persistent high/hard 
state (unknown in Galactic BHBs), which can extend all the 
way to the Eddington luminosity and smoothly morph into the 
very high state or other super-Eddington accretion states. 
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Figure 2: A typical ULX spectral shape (grey) in 0.3-10 keV, with a soft ex- 
cess below 2 keV and a hard curvature above 2 keV. The soft excess can be 
modeled by a cool thermal component, over the power-law extension of the 
hard component. A slim disk model (p-free) or a warm, thick Comptonization 
model can adequately fit the hard curvature; in the plot, a p-free model with 
kT m = 2.5 keV and p = 0.5 has an almost identical spectral shape to a compTT 
model with kT e = 2.04 keV and r = 6.39, which is shifted vertically for clarity. 
For comparison, two diskbb models are plotted with kT m = 1.0 and 2.5 keV, 
respectively. They are more curved and cannot fit the energy spectrum over the 
0.3-10 keV band. 



At the opposite end of the broad photon-index distribution, 
several ULXs with T > 2.5 could be classified in the steep 
power-law state of Galactic BHBs, typically f o und at near- 



Eddington luminosities (IFeng & Kaaretl l2005t IWinter et al 



2006; Soria 2007). The fact that there are no clear gaps, but in- 
stead a continuous distribution in photon index and an overlap- 
ping luminosity distribution between hard and soft ULXs may 
be an indication that the two regimes are smoothly connected. 

In conclusion, with longer X-ray observations and increasing 
photon statistics, we have reached a stage where we can identify 
and model a number of subtle spectral features on top of the 
simple power-law for many ULXs. In the rest of this section, 
we discuss a variety of spectral models that account for such 
features, and their physical implications. 

4.2. Soft excesses: cool disk emission or massive outflow? 

In many ULXs, excessive residuals at low energies (< 2 keV) 
would be seen if one makes a power-law fit at high energies; 
an additional disk component could adequately fit the excesses 
and the inferred disk inner temperature is arou nd 0.1-0.4 keV. 
This c ool, soft excess was first identified bv | Fabian & W ard 
d 1993b in NGC 5408 X-l; iKaaret et al.1 d2003h showed the 
emission was point like and applied a disk plus power-law 
model. Simil ar phenomena we re then found in NGC 1313 
X-l an d X-2 dMiller et al.ll2003|). and later on in many other 



ULXs dMiUer et al 120041: iFeng & Kaaretl l2005t IStobbart et al" 
20061: IWinteretal.1 12006). Such a fit is analogous to the em- 



pirical model for Galactic BHBs except that the disk is cooler 
and more luminous in ULXs than in sMBHs at their high 
state (Gierlihski & Done 2004). For a standard accretion disk 



(Sha kura & SunvaevNl973l) extending to the last stable orbit 
around a BH, the inner disk temperature T\ B is a function of 
BH mass M and total luminosity L: 



1/4 



I- 



M 



\L E dd) \1OM 



-1/4 



keV, 



(3) 

where a depends on the BH spin (a = 1 for a Schwarzschild 
BH, a — 1/6 for an extreme Kerr BH), £ is a corrective factor 
that takes into account the no-torque conditio n at the innermost 



orbit , and k is the spectral hardening factor (Makishima et al 



2000). Given the same spinning and correction factors, the BH 
mass is scaled with the disk inner temperature and luminosity 
as a function of 

M oc T in 2 L 1/2 (4) 

Thus, the low-temperature, high-luminosity characters of the 
assumed disk suggest that the emission arises from accretion 
onto IMBHs. For comparison, the typical disk inner tempera- 
ture is about 1 keV at a luminosity of 10 38 erg s _1 for sMBHs 
about 10 M . The characteristic temperature for ULXs is usu- 
ally around 0.1-0.4 keV and the luminosity is some 10 39 - 
10 40 erg s _1 , suggestive of a BH mass of 10 3 - 10 4 M Q , 

However, the interpretation of soft excesses as emis- 
sion from cool a c cretio n disk has confronted challenges. 
Goncalves & S oria 2006) argued that soft excesses could also 



be soft deficits depending on the e nergy band where the power- 
law is fitted. Berg hea et al.1 (12008b found similar soft excesses 
in both ULXs and less luminous X-ray sources, implying that 
they are not a unique signature of IMBHs in a statistical sense. 
It is suggested that supercritical accretion may produce mas- 
sive outflows, which are optically thick with soft, thermal emis- 
sion at similar temperature , being ano t her possible origin for 
the detected soft exc esses {King 2004; Begelman et al. 2006; 
I Poutenen etail l2007h Other authors dRobertsI 120071: ISoria 
l2007t iGladstone et al.1 120091: IWeng & Zhang||201 lb have sug- 
gested that the soft excess does come from the disk surface, but 
much further away than the innermost stable orbit, where the 
disk may be obscured or replaced by an outflow or a scattering 
corona (Section [4. 3l l. In that case, the characteristic radius of 
the thermally-emitting disk region cannot be used to infer the 
BH mass. 

Repeated observati ons of an individual ULX may h elp test 
the cool disk model ( Soria 2007 ; Feng & Kaaretll2007bb . If the 
disk touches the innermost stable orbit around the BH, one ex- 
pects to see the disk luminosity vary to the 4th power of the in- 
ner temperature, L jisk K , as well known in the thermal state 
of Galactic BHBs (Gierlihski & Donel20 04). If, instead, the ob- 
served inner disk radius is the boundary between an outer stan- 
dard disk and an inner corona or outflow, we expect that such 
transition will move further and further out at higher accretion 
rates whe n super-Eddington. In t he simplest outflow sc enario, 
Ri n ~ >h (Begelman et al. 2006; Pout anen et al. 12007b . This 
implies dSoriall2007b Ldisk ~ constant for a viscously heated 
outer disk and Ldisk K ^ in 4 for a strongly X-ray-irradiated outer 
disk (as is more likely the case in ULXs). Most other plau- 
sible relations between R m and tin predict negative correlations 
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between the disk luminosity and temp erature, with slopes be- 
tween and -4 jPoutanen et alj |2007). A negative correlation 



between the luminosity and temperatur e is clearly detected for 
at lea st NGC 1313 X-2 and IC 342 X-l dFeng & Kaaret|l2007bl 
2009h . From a sample of 8 ULXs, iKaiava & Poutanenl T2009) 
found that the luminosity of the soft excesses scales with the 
temperature as L so f t oc r~ 3 5 , roughly consistent with the model 
prediction that the soft emission arises from a photosphere due 
to strong outflows from a supercritical disk a nd its tempera 



hire is determined by the spherization radius ( Poutanen et al 



2007; Soria 2007). On the contrary, a positive correlation is 
seen from N GC 5204 X-l; howeve r, it has a shallowe r slope 



McClintock et al. (2009) 



L cc ;r.2-i±o.5 (|Feng & Kaaretl 120091) 
suggested, based on a disk atmosphere model, that the harden- 
ing correction factor becomes larger at increasing luminosity, 
and explains the apparent shallower slope in Galactic BHs such 
as XTE J1550-564 and H1743-322. Applying such a varying 
hardening correction, one can recover the L oc T A relation for 
NGC 5204 X-l, which would favo r the IMBH scenario How- 
ever, this result is not robust yet: IKaiava & Poutanen! ([2009) 
found that the positive correlation disappears or even becomes 
negative if one chooses a different lower energy bound for the 
spectral fitting. Also, the disk contributes only < 25% of the 
0.3-10 keV luminosity, inconsistent with the thermal dominant 
state. 

Instead, if ULXs were IMBHs, we would expect to see many 
of them in the canonical thermal state dominated by the soft 
disk component. In summary, using the disk parameters in- 
ferred from X-ray spectra to determine the BH mass is based 
on assumptions that have proved either incorrect or problem- 
atic for most sources. 



4.3. High energy curvature: hot standard disk, slim disk, or 
warm corona ? 



It w as noted by IStobbart et al. I d2006h and iGladston e et al 
d2009l) that high energy curvature/steepening is a com- 
mon feature in high-quality ULX spectra. Generally, the 
break/steepening could be fitted by a Comptonization model 
with relatively low energy electr ons and high optical depth 
dMakishimal 120071; iRobertsI 120071) . A few ULX spectra are 
mildly convex, formally consistent with a standard disk model. 
We thus discuss the physical implications of these models. 

4.3.1. Hot standard disk and thermal state 

With ASCA observations, a hot multicolor accretion disk 
model was found to fit the spectra of a number of ULXs {e.g. 
Kub ota et alJ 12001). For those with relatively low tempera- 



tures (kT m < 1 keV), t hey could be expl ained by accretion 
onto sMBHs or MsBHs d Winter et alJl2006h . The high tem- 
perature and luminosity may be explained by large BH spin 
(Ma kishima et al J 120001) , but this does not work for extreme 
cases where relatively small BH mass is required, which in 
turn is inconsistent with the observed high luminosity. For 
XMM-Newton spectra, a single hot disk with absorption is usu- 
ally unable to fit the data; an additional co mponent is required 
and dominates in the low energy band dRoberts et alJ 12005 



Feng & Kaaretll2005t IStobbart et al.ll2006t IWinter et ailEoO 
and th e disk component is not the dominant part dStobbart et al 
2006). The presence of a soft component in addition to the hot 
disk is unusual in BHBs and does not have a simple physical 
interpretation. In general, the hot disk model is problematic. 

The standard disk emission is best known in the thermal state 
in BHBs, in which the X-ray spectrum is dominated by emis- 
sion from a standard accretio n disk; plus, the X-ray variabilit y 
is low, with no or weak QPOs dRemillard & McClintockl2 006). 
So far, the thermal state may have been identified only in one 
case, M82 X-l, constrained using both spectral and timing in- 
formation. With si multaneous Chandra an d XMM-Newton ob- 
servations of M82, Fen g & Kaaret (2010) found surprisingly 
that the previously known QPOs in the source disappeared. 
The light curve was no longer highly variable and the power 
spectrum was consistent with that of white noise. The energy 
spectrum also changed dramatically from a straight power-law 
to a disk-like spectrum. Fitting with a standard accretion disk 
model, the inferred luminosity and temperature from three re- 
peated observations are consistent with a 4th power relation. 
All results are well consistent with that expected for the ther- 
mal state. The monitoring data from RXTE indicate that these 
Chandra and XMM-Newton observations were made during the 
source outbursts, suggesting that M82 X-l usually stays in the 
hard state and could transition to the thermal state during out- 
bursts. Fitting with a relativistic disk model, the BH is esti- 
mated to be a fast spinning IMBH of 200-8 00 M Q 



There are a few more sources dlsobe et al. 2008 



Kajava & Poutanen 2009; Jin etaDEoiol) which have shown 
spectral variability possibly consistent with a L oc T 4 relation 
if fitted with a hot standard disk model. Another group of five 
nearby ULXs consistent with the thermal state sp ectrum, based 
on X- ray spectral modeling, was discussed by iHui & Krolik 
(2008). These sources all lie at the lower end of the ULX 
luminosity regime, and could be explained by Kerr MsBHs 
in the thermal state. They do not stand for the majority of 
ULXs, in particular those with a typical luminosity around 
10 40 erg s _1 . Also, lack of fast timing information hampers a 
firm identification. 

Therefore, it is gen erally accepted th at the canonical thermal 
state is rare in ULXs dSoria et al J l2009h . 



4.3.2. Slim disk and p-free model 

The slim disk can radiate near or moderately above 
the Eddington limit estimated from spherical accretion 
dAbramowicz et al" 1988b . This super- Ed dington disk i s thus 



a good candidate for ULX emission dWatarai et al.l 12001 



Ebis awa et alJl2003l) . When advective energy transport dom- 



inates over radiative cooling in an accretion disk, the ra- 
dial temperature profile changes from T oc j?~ 3 ^ 4 for a stan - 
dard disk to T cc /r 1/2 for a slim disk dWatarai et al.ll2000h . 
An extended multicolor disk model called the p-free model 
was introduced, with the radial dependence p on temperature 
(T oc R~ p ) being an additional free parame ter dMineshige et al] 
1994 ), and used to expla in ULX emission dWatarai et al.ll2.Q0lh . 
Vierdayanti et al. I d2006l) investigated the energy spectra of four 
ULXs, NGC 5204 X-l, NGC 4559 X-7/X-10, and NGC 1313 
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X-2, and found that they can be well fitted by the p-free model. 
Interestingly, the best-fit p value was found to be around 0.5, 
consistent with that for a slim disk. In the luminosity vs. disk 
temperature diagram, the four ULXs and Galactic BHBs ap- 
pear to be in line with the same L oc r 4 relation, suggesting 
that ULXs may be BHs of similar or slightly higher masses 
(sMBHs or MsBHs) but in the slim disk state. Another suc- 
cessful test of the p-free or slim disk model was for M33 X- 
8, which is marginally classified as a ULX but a good candi- 
date for deep investigation owing to its relatively sma ll distance 



(e.g., iFoschini et al.l l2004. 20061) . IWeng et alJ (120091) first iden 



tified that an extra power-law component is required in addition 
to the p-free disk to fit the source spectrum. The best-fit p value 
(0.57 + 0.03) suggests notable advection in the disk, close to the 
slim disk nature rather than a standard one. These indicate that 
the slim disk state may be a natural extension of the canonical 
thermal state at supercritical accretion. This is consistent with 
the radiation-magnetohydrodynamic simulations that the hard, 
thermal and slim states could be produced under a single model 
with different accretion rates dOhsuga et al.ll2009h . 

However, the p-free model encounters difficulties in both 
observation and theory. In some sources, the disk inner tem- 
perature was found too high t o be physically reason able when 
fitted with the p-free model (Gla dstone et al. 2009). If pho- 
ton trapping is taken into account, the radial temperature pro- 
file is expected to b e significantly flatte r than p = 0.5 within 
the trapping radius dOhsuga et al . 2002, 2005), in contrast to 
the spectrum predi cted without considering photon trapping 
dWatarai et al. 2000). Also, numerical simulations suggest that 
dense, wa rm outflows are always a ssociated with supercritical 
accretion ( Ohsuga et al J2005ll2009l) . which would Comptonize 
and obscure the emission from the underlying disk, see details 
in the following subsection. 

4.3.3. Warm corona and ultraluminous state 

In Galactic BHBs, high energy spectral curvature is also de- 
tected, but usually occurs at tens to hundreds of keV. It is in- 
terpreted as Comptonization in a thin corona with temperatures 
typical of 10 2 keV. However, the spectral curvature in ULXs 
suggests an electron temperature of a few keV and a scatter- 
ing optical depth significantly higher than unity, indicative of 
a warm and thick corona. Due to the presence of such a thick 
corona, the emission from the disk would be distorted and the 
previous interp retation, a simp le summation of a cool disk and 
a warm corona ( Stobbart et al.ll2006l) . needs be modified conse- 
q uently. 

Gladstone et al. d2009l) revisited best quality XMM-Newton 



data and concluded that two common features simultaneously 
exist in the energy spectrum of ULXs: soft excesses below 2 
keV and high energy curvature above 2 keV. This is not fre- 
quently seen in the spectra of Galactic BHBs when they are 
thought to accrete at sub-Eddington level. Thus, they suggested 
that ULXs represent a new state distinct from the four well- 
known states (quie scent, hard, thermal, and steep power-law) 
in Galactic BHBs dRemillard & McChntockll2006l). and ca ll it 
ultraluminous state. Please refer to Rob erts et all d2010l) for 
a more detailed review of this state. The Galactic BHB that 



comes closest to this state is XTE J 1550-564 at the peak of 
its 1998 September outburst . A coupled disk-coro na model 
dSvensson & Zdziarskil ll994; Don e & Kubota l2006h was pro- 
posed to explain such state, in which the disk and corona are 
not independent but share the total gravitational energy release. 
The un-scattered outer disk is responsible for the soft excesses, 
while the Comptonized inner disk produces the hard, curved 
tail. The BH mass, estimated from the unmasked temperature 
and radius of the underlying disk, lies in the regime of sMBHs 
or MsBHs. 

Numerical simulations of supercritical accretion indicates 
that radiation pressure w ould power strong outflow s above the 
disk in the inner region dOhsuga et alJl2005l 120091) . The wind 
is optically thick and warm (a few keV), and would strongly 
Comptonize the disk emission. At high accretion rate, the disk 
becomes thick and would collimate the radiation along the sym- 
metric axis causing geometric, mild beaming with an appar- 
ent isotropic luminosity over the Eddington limit by a factor 
up to ~10. This picture is consistent with the proposed ultra- 
luminous state. Observationally, f rom the spectral evolution of 
M33 X-8. lMiddleton et all d2011bl) suggested that the outflow is 
launched at near-Eddington luminosity, and the radius at which 
it is launched moves further out with increasing luminosity. 

4.3.4. Bulk motion Comptonization 

If the Comptonizing region is the fast outflow, rather 
than a static corona, it is possible that bulk motion 
( Titarchuk & Zanniasl 1998 ) can be as or more important than 
thermal Comp tonization, depending on the characteristics of 
the outflow ( Kawashima et alJ l2Q09h . Observationally, ULX 
spectra that can be formally fitted with thermal Comptonization 
models can be equa lly well fitted with bulk motion models {e.g. 
Rober ts et al 120051) . The only difference is in the interpretation 
of the c haracteristic energy of the upscattering electrons. As 
noted bv lSorial d201 lb . simple physical considerations, compar- 
isons with other astrophysical classes of outflows, and numeri- 
cal simulations all suggest that super-Eddington ULX outflows 
have a characteristic speed ~ 0.1c ~ escape velocity from the 
launching radius, corresponding to a bulk kinetic energy ~ 5 
keV for the outflowing electrons. This is similar to the charac- 
teristic energy of the spectral breaks. 

4.4. Supersoft ULXs 

Most ULXs and BHBs have considerable radiation output 
above 2 keV, while supersoft sources (SSSs) have most of their 
energy released below 2 keV. Canonical SSSs are speculated to 
be pow ered by steady nuclea r burning on the surface of white 
dwarfs (IDi Stefano et alJl2010h . A number of supersoft ULXs 
have been discovered, whose spectrum is dominated by a cool 
blackbody component of tens to a hundred eV, being a distinct 
cla ss from the majori t y of ULXs. These include M101 ULX 



1 dPence et all 1200 it Kong & Pi Stefano) 120 05: Mu kai et al 



2005 ), Antennae X - 13 dFabbiano et alJ (20031). M81 ULS1 



2011 



(ISwartz et al.ll2002t iLiul 120 08). NGC 4631 XI (ICarpano et al 
20071: ISoria & Ghoshl 12009k and NGC 247 ULX djineta" 



). At long timescales, most of these sources have shown 
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Figure 3: A toy model of ULX. The source is powered by supercritical ac- 
cretion onto a black hole. Radiation pressure propels strong outflows around 
the inner region of the accretion flow, where a slim disk may exist in connec- 
tion with a standard disk far outside. The outflow is warm (a few keV) and 
optical-thick and may scatter off and redirect photons, framing a funnel with 
geometric beaming; it could also up-Comptonize soft photons; the inner disk 
and/or the outflow are responsible for the observed hard emission component 
above 2 keV with a curvature. The soft excess may arise from the outer disk 
or the optical-thick outflow. If viewing along the disk plane, one may only see 
the soft emission of the outflow and this could be one scenario for the supersoft 
ULXs. When the outflow hits the interstellar medium, it may produce shocks 
seen as optical nebulae. The picture is not to scale. 



dramatic variability with a luminosity change by a factor up 
to 10 3 but with consistently low temperatures. A wea k power- 
law ta il is significantly detected in at least one source djin et al.l 
201 1 ). but none of the supersoft ULXs has become a power-law 
dominated ULX or vice versa. Fast timing noise w as detected 
in two of them, M101 U LX-1 dMukai et al.l l2005h and NGC 
247 ULX djin et alJl20lH) . down to a timescale around 100 s, 
indicating that they are not AGN. The nature of these peculiar 
sources is very uncertain. They are unlikely powered by white 



4.5. State transitions of ULXs 

The study of long term spectral variability of ULXs mainly 
relies on pointed observations with Chandra or XMM-Newton, 
or quasi-regular snapshots with the X-Ray Telescope on board 
Swift. In addition, RXTE has been used to monitor the brightest 
(in terms of observed flux) ULX in M82. The most interest- 
ing science with these observations, in addition to the search 
for orbital periods, is perhaps to find state transitions known in 
Galactic BHBs. 

Canonical transitions between the low/hard and high/soft 
states were reported in IC 34 2 X-l and X-2 with ASCA ob- 
servations ( Kubota et alj|200l|) and also in Ho lmberg IX X-l 



with multiple missions (ILa Parola e t al. 2001). However, al 



most no such canonical state transitions have been reported in 
the Chandra and XMM-Newton era. In contrast, transitions be- 
tween soft and hard sta tes were found but not in a way simi- 
lar to Galactic sources (|Dewangan et al.l 120041 iFeng & Kaaretl 
2006bl 120091: lDewanganetalJl2010l) . A. major difference be- 
tween ULXs and Galactic BHBs is that most ULXs are persis- 
tent sources and have been bright since their discovery decades 
ago, while the latter are usually transients especially accreting 

from a low mass co mpanion with Roche lo b e overflow. 

Swift monitoring dKaaret & Fend 12 009: Strohmaver 2009; 
Kong et alJl2O10t IVierdavanti et alJbOlOt iGrise et alj|2010l) of 



bright ULXs have revealed flaring activity with a luminosity 
change by a factor of < 10 on tim e scales of days to months. In 
particular, iKaaret & Fengl |2009) found Ho lmberg IX X-l re- 



mains in the hard state as its flux varies, and Grise et al. (2010) 



found the X-ray spectral state is not correlated with the lumi- 
nosity in H olmberg II X- 1 . Also, as mentioned previously in 
Section 1431 iFeng & Kaarel d201dh found M82 X-l may have 
changed from the hard state to thermal state during outbursts. 

There is one exceptional case. ES0243-49 HLX-1 shows 
a low/hard to high/soft jump reminiscent of the c anonical state 
transition dGodet et alJl2009l: IServillat et al.ll201 lb . It also dis- 
plays outbursts in a "fast rise exponential d ecay" pattern, ver y 
similar to those seen in Galactic BHBs dLasota et al. 201 1 ). 
However, the authors argued that the observed variability is un- 
likely triggered by disk instability due to inconsistent timescale. 



dwarfs, because they are too luminous (a few times 10 y erg s~') 5. Fast X-ray variability and mass scaling 
and some of them are too hot (> 100 eV). They are also un- 
likely to be due to cool disks around IMBHs, because such disks 
cannot have a huge variability with a constant temperature. One 
possible scenario is that their emission is produced by massive 
outflows from supe rcritical accreting sources viewed at a high 
inclination angle ( Ohsu gaet al.ll2005 : Poutanen et al. 2007 ). A 
unification model (supercritical accretion with thick outflows, 
see Figure 0) may be able to explain both normal, broad-band 
ULXs and less luminous, supersoft ULXs with different view- 
ing angles. 

So far, we have briefly reviewed a variety of scenarios that 
can account for the observed ULX spectra (see also Figure |2}. 
In Figure[3] we sketch the physical ingredients behind the most 
common combination of such models, which may be success- 
fully used to explain the multiband spectrum of most ULXs. 



5.1. Observational characteristics 

The power spectrum of Galactic BHBs can be decomposed 
into broad components (power-law red noise, band-limited 
noise, flat-top or zero-centered broad Lorentzian, etc.) and 
narrow components (QPOs). Although we still know little 
about the physical origin of these temporal features, their phe- 
nomenology h as proved use ful in defining and identifying ac- 
cretion states dBellonil l2010) and possibly in determining BH 
masse s via model independent calibrations (e.g. lMcHardv et all 
l2006h . 

An important model-independent finding of iHeil et al 



(2009), based on the study of sample of 16 bright ULXs, is 
that those sources can be divided into two groups defined by 
their short-term X-ray variability, even though they have similar 
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energy spectra and luminosities. One group (including among 
others NGC 5204 X-l) shows weak or absent variability on 
timescales at least as large as 100s. The other group (includ- 
ing among others NGC 5408 X-l) has similar variability levels 
and power spectra as luminous Galactic BHBs and AGN in the 
observed frequency bandpass (10~ 3 -1 Hz). It is still unclear 
what suppresses short-term variability in one class but not the 
other. In Galactic BHBs, high variability is often associated to 
non-thermal emission wit h steady jets; l ow variability to ther- 
mal emission with no jets dBellonfcOlOh . 

In the following, we select four ULXs in which interesting 
structures have been found in their power spectra, shown in Fig- 
ure m 

M82 X-l, also known as M82 X41.4+60 — 



Strohmaver & Mushotzkvl Q2003) first discovered low fre- 
quency QPOs around 54 mHz from the M82 galaxy with 
both XMM-Newton and RXTE observations. By fitting the 
XMM-Newton images to its poin t spread function with k nown 
source positions from Chandra, IFeng & Ka aret (2007a) iden- 
tified that M82 X-l is the source that produced the QPOs, and 
interestingly, it was not the brightest X-ray source in M82 at the 
time of QPO discovery. A following XMM-Newton observation 
revealed that the QPOs have changed the central frequency 
from 54 mHz to 114 mHz, above a flat-top component with 



a low frequency break at 34 mHz ( Dewanga n et al. 120061: 



Mucciarelli et a i1 l2006h . Mucciarelli et al. (2006) also found 



that the QPOs changed frequency from 107 mHz to 120 mHz 
during the second observation, and a plausible harmonic 
ratio of 1:2:3 for QPO frequencies with RXTE data. These 
QPOs, however, were undetected with three XMM-Newton 
observatio ns of similar depth, suggestive of a spectral state 
transition (IFeng & KaarefeOlOl) . 

NGC 5408 X-l — With a short XMM-Newton observation, 
Soria & Motchl ( 2004 ) found the source is hi ghly variable, and 



its po wer spectrum has a break at 2.5 mHz. IStrohmaver et al 



(2007) reported the identification of narrow QPOs around 
20 mHz from the source, and confirmed the presence of the fre- 
quency break. The QPOs are narrow and may show a pair with 
a 4:3 frequency ratio. A longer XMM-Newton observation later 
detected the QPOs at about 10 mHz and the energy spectral pa- 
ramet ers of the source also changed dStrohmaver & Mushotzkv 
2009). The QPO frequency was found to scale with the disk 
flux and also the spectral index of the power-law component. 
These correlations are analogous to that seen in Galactic BHBs, 
suggesting that this ULX is in an emission state similar to sub- 
Eddington sources. The QPO properties (narrow, strong, with 
a flat-top continuum) and the temporal-spectral evolution pat- 
terns suggest that the QPOs are of type C and the compact ob- 
ject in NGC 5408 X-l is an IMBH. However, as the QPO fre- 
quency is not proportional to the break frequency of the contin- 
uum as expected, the identificatio n as of type C was ques tioned 
and so was the IMBH scenario (Mid dleton et alj|201 lal) : they 
proposed that the source is super-Eddington and the turbulence 
in the Comptonization region gives rise to the observed vari- 
ability. Oth er authors have also rejected the IMBH argument 
dSoriall2007h . pointing out that neither lower QPO frequencies 
nor larger, cooler disks can be used as indicators of higher BH 



masses in ULXs; both findings may simply suggest that the vis- 
ible inner radius of the disk has receded much further than the 
innermost stable orbit, and has been covered or replaced by a 
thick Comptonizing region o r outflow in the inner region. Be- 
sides, Heil & Vaughanl l 20101) also found that the source shows 
a linear correlation between the rms amplitude and flux, similar 
to that seen in bright BHBs and AGN. 

M82 X42.3+59 — This is a transient ULX in M82 



(Mat sumoto et alj |2001). It could not be detected by Chandra 
in its low state, but sometimes, was even bri ghter than M82 
X-l, being the most luminous source in M82 (Feng & K aaret 
2007a). With Chandra observations, and also confirmed by a si- 
multaneous XMM-Newton observation once, IFeng etalJd2010l) 
discovered low frequency QPOs around 3-4 mHz in its power 
spectrum. These QPOs are wide and strong, appear above white 
noise, and can only be detected when the source is brighter than 
10 40 erg s _1 . This is consistent with a Type A or B classifica- 
tion, but is inconsistent with the properties of type C QPOs. 

NGC 6946 X-l — This source has the highest short-term 
variability found in a ULX, with a fractional rms am plitude of 
60% i ntegrated in the frequency range of 1-100 mHz dRao etalJ 
120101) . Its power spectrum shows a flat-top continuum that 
breaks at about 3 mHz with possible quasi-periodic oscillations 
(QPOs) near 8.5 mHz. 

5.2. Implications on the mass 

The characteristic frequency alone usually is not a direct in- 
dicator of the BH mass, as these noise features can vary dra- 
matically in a single source. However, correlations between 
the frequency and spectral parameters have been found, and the 
correlation patterns are thought to scale with the BH mass. 

T he most common QPOs seen in BHBs are of type C 
(cf. ICasella et all 120051) . The frequency of type C QPOs 
are found to vary over a wide range. For example, the 
lowest d etected frequencies in GRS 1915+105 are around 
1 19971) . even lower than the frequen- 



1 mHz ( Morg an et al. 



cies of all QPOs detected in ULXs. A positive correlation 
between the QPO frequency and disk flux has been found 
for type C QPOs in XTE J 1550-5 64, GRO J 1655 -40 and 
17 43-322 dSobczaketalj EoOol: " 
t al.l l2009h . The 



H 



Remillard & McClintock 



2006; lMcClintock et al.ll2009l) . The QPO frequency also scales 
with the photon in dex of the power-law component in the en - 
ergy spectrum (e.g. Sobczaketal. 2000; Vi gnarca et al. 2003). 
Titarclmk & FioritoT i 2004h suggested that the QPO frequency 



scales inversely with the BH mass at a fixed power-law pho- 
ton i ndex, and thus one can infer the mass given spectral in- 
dex dShaposhnikov & Titarchuk 2009). On the contrary, Type 
A and B QPOs show in a relatively narrow frequency range. 

The broad continuum seems to be a more reliable deter- 
minant of BH mass, as it is detected in both BHBs and 
AGN and the scaling has been tested over a wide frequency 
and mass range. sMBHs and AGN appear to populate a 
"variability plane" defined by the BH mass, accretion rate, 
and v/, which is the frequency of the lower high-frequency 
Lorentzian component in the power spectrum (McHardv et al] 
2006; iKbrding et al 1 l2007h . Limited by sensitivity, v/ cannot be 



detected in ULXs. Based on a correlation between v/ and the 
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Figure 4: Power spectr a of four ULXs. M82 X-l: reproduced using data from 
iFeng & Kaarel <2007al) with XMM-Newton observations in 2001 and 2004, re- 
spectively; A and B refer to the photon extraction regions. NGC 5408 X-l: the 
two spectra are respectively from XMM-New ton observations in 2006 and 2008 
(Figure 5 in Stro hmaver & Mu shotzkv 2009). NGC 6946 X-l: reproduced us- 
ing data from Rao et al. 1 2010). M82 X42.3+59: four spectra are from Chandra 
(b-d) and XMM-Newton (a) observations (Figure 1 in Fens et al. 2 oTp|) . 



frequency of type C QPOs, ICasella et al.l ([2008 ) extended the 
variability plane to type C QPOs that are visible in ULXs. The 
BH mass can be inferred from the X-ray luminosity (here we 
adopt the radiatively efficient accretion case) and type C QPO 
frequency as 

log(M/M ) > 0.5 log(L x /0.1c 2 ) - 0.51 log v QP o - 7.9 , (5) 

where the inequality is due to the unknown bolometric correc- 
tion (Lx ;$ ^boi)- Although the "variability plane" is tested 
for a wide range of BH masses, the correlation between v/ 
and vqpo has not been tested in ULXs and even the identifi- 
cation as of type C is q uestioned at least for NGC 5408 X-l 
( Middleto netai1l201 lab . 



Attempts to scale BH mass with the varia bility amplitud e 



(rms) have made great success in AGN (IZhou et 



amplitude 

IDEoio). 



However, it is found that ULXs are significantly less variable 
than AGN, and a simple ext rapolation could not be obtained 
dGonzalez Martin et al.ll201 lb . 

Based on the correlations mentioned above, the BH masses 
inferred from the characteristic frequencies are all about 10 3 - 
10 4 M , indicating that ULXs contain IMBHs. However, appli- 
cation of these mass scaling methods requires strong caveats. 
There are at least three problems. First, we are not sure whether 
the timing features detected in ULXs are the same type found 



to scale with BH masses in Galactic BHBs; for example, we 
are not sure whether the Q POs in NGC 5408 X-l are of Type 
C (IMiddleton et alj|201 lab . For other types of timing features 
{e.g., type A and B QPOs), it is still not k nown whether thei r 
frequency is proportional to the BH mass (ICasella et al ] l2005h . 
Second, frequency-luminosity-mass relations have been cali- 
brated and tested for Galactic BHBs and AGN in sub-Eddington 
accretion, that is when the inner edge of the disk is constant or 
moves towards the innermost stable orbit as the accretion rate 
increases (thus, increasing QPO frequencies). There is no em- 
pirical evidence that the scaling goes the same way in ULXs, 
which may be super-Eddington systems where the visible inner 
edge of the disk recedes from the innermost stable orbit for in- 
creasing accretion rates (thus, reducing QPO frequencies). The 
third problem is that we do not have enough signal-to-noise ra- 
tio to detect high-frequency QPOs, which severely limits possi- 
bility of direct comparisons with Galactic BHBs. 

6. Optical counterparts and stellar environment 

In this section the "optical counterpart" refers to the point- 
like optical source that is spatially associated with the ULX; 
sources with resolved structure, such as optical and radio nebu- 
lae, will be discussed in the next section. The optical emission 
from the counterpart could arise from the donor star or outer 
accretion disk, or both. It provides interesting information re- 
garding the binary evolution history, nature of the donor star, 
disk geometry, and mode of mass transfer, and would also help 
constrain the BH mass. 

Several groups have modeled X-ray irradiation in order 
to distinguish between the contributions of the outer disk 
and th e donor star, and constrain the intrinsic stellar prop- 
erties. ICopperwheat et al.l (120051 120071) applied their irradia- 



tion model to the photometric data for seven ULXs; they sug- 
gested that donor stars are consistent with spectral type B or 
later, and the BH masses are in the range of ~ 10-100M o . 
Patru no & Zampier i (2008, 2010) computed evolution tracks of 
ULX counterparts on the color-magnitude diagram, assuming 
a massive companion orbiting an sMBH or a MsBH; they ap- 
plied their models to NGC 1313 X-2, suggesting a BH mass 
~ 50-IOOMq. In addition to reproducing the observed colors 
and luminosity (after accounting for irradiation and disk contri- 
bution), a donor star in a ULX must be able to supply a mass 
accretion rate > 1O~ 6 M yr~', which requires Roche lobe over- 
flow, and to do it over a sufficiently long period of time to ex- 
plain the observed number of ULXs per galaxy. Binary pop- 
ulation synthesis m odels have been developed and used to ad- 
dress this problem dPod siadlows ki et al.ll2003l: iRappaport et al. 
l2005t iMadhusudhan et al.ll2006l 120081) . Modeling the binary 
evolution requires calculations of the mass transfer rate, the 
corresponding change in orbital separation, and the subsequent 
fate of the transferred mass. This in turn depends critically on 
whether we impose a strict Eddington limit on the mass trans- 
fer rate that can be accepted by the BH (with the rest being 
expelled in a wind), or instead we assume that the BH can ac- 
crete at higher rates. Thus, for the same kind of donor star, 
the outcome in terms of X-ray luminosity and lifetime is very 
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different depending on whether the primary is an Eddington- 
limited stellar BH, a non-Eddington-limited stellar BH, or an 
IMBH. 



Madhusudhan et al 



20081) considered three character- 
istic families of binary evolution tracks, with different BH mass 
distributions: 6-15M , 6-24M , and a single-mass population 
of 1OOOM BHs. They built probability diagrams in different 
slices of parameter space (optical color-magnitude diagram; or- 
bital period versus potential X-ray luminosity; X-ray luminos- 
ity versus age), representing the total time spent by all the sys- 
tems in those regions of parameter space. Their main conclu- 
sion is that the most probable ULX system parameters corre- 
spond to high-mass donors (> 25M ), with orbital period be- 
tween 1 and 10 days and characteristic age ~ 10 7 yr. In terms of 
BH mass, they found that both the 1OOOM BH population, and 
that with non-Eddington-limited BHs up 25M were consistent 
with the optical and X-ray observations and with the observed 
number of ULXs per galaxy; the 25-lOOM BH mass range has 
not been tested yet. 

Thanks to the sub-arcsecond resolution of Chandra and HST, 
the optical counterparts of more than a dozen ULXs have been 
identified. These are done mainly by aligning Chandra and 
HST images using objects with emission on both, to reduce the 
Chandra X-ra y position error from 0.6" to about 0.2"-0.4". 



Tao et al 



201 11) collected all available data in the HST archive 
to investigate the optical properties of 13 ULXs, which can be 
summarized as follows. 



a) The X-ray to optical flux ratios and optical colors are consis- 
tent with that of LMXBs, but unlike HMXBs or AGN. This 
indicates that the optical emission from ULXs is mainly due 
to reprocessing of X-rays on the outer accretion disk. As a 
consequence, one cannot infer the spectral type of the com- 
panion star using optical colors directly. 

b) The broad band optical spectrum can be generally fitted with 
a power-law model (F v oc y"). The distribution of the spec- 
tral index a peaks between 1 and 2, which is also consistent 
with disk irradiation models. 

c) The counterpart has unique optical properties within its 
neighborhood. Nearby sources usually show a blackbody 
spectrum instead of a power-law spectrum. 

d) The optical emitting region is estimated to be large, of 
the order of 10 12 cm, indicative of long binary periods 
(> 10 days) for sMBH or MsBH primaries, or shorter peri- 
ods for IMBHs. 

e) Sources with multiple observations usually show strong, 
random variability on both magnitudes and colors. 

There are a few outliers showing peculiar properties. The 
optical spectra of NGC 2403 X-l and M83 IXO 82 have a slope 
consistent with that of the standard disk emission, F v oc y 1/3 . 
M101 ULX-1 and M81 ULS1 are more similar to HMXBs than 
LMXBs in terms of X-ray to optical flux ratios. IC 342 X-l 
may show a blackbody-like spectrum inste ad of a power-law. 
More details can be found in Tao et alJ ( 2011 ). 

There have been a lot of observational efforts to determine 
the nature of the companion stars in nearby ULXs. So far, NGC 
7793 PI 3 is the only ULX in which emission from the compan- 
ion star is directly seen beyond reasonable doubt dMotch et al 



201 II) : detection of narrow emission lines identifies the com- 
panion star as a late-B type supergiant with a mass « 2OM . 



Orbital periods may have been detected in a few sources. 
Assuming Roche-lobe overflow, the companion density can be 
simply determined from the binary period, p — 1 15P hr 2 g cm" 3 
(Fra nk et alJ |2002). M82 X-l is the strongest case of a 
positive detection; a monitoring progr am with RXTE r e- 
vealed a steady perio d of about 62 day dKaaret et alj|2006albl : 



Kaaret & Feng 2007), corresponding to a companion density of 



5x10 5 pcm 3 Jstrohmavei ( 20091) discovered a 1 15 day X-ray 



period from NGC 5408 X-l with Swift monitoring, and inferred 
a similar companion density of 1.5 x 10" 5 g cm" 3 ; however, the 
interpreta tion of that modu lation as a binary period has been 
disputed dFoster et al. 2010h . These claimed long periods im- 
ply an evolved donor star. Possible detections include a 6.1 day 



perio d from NGC 1313 X-2 with HST observations dLiu et al 



120091) : howe ver, low significance of the period was claimed b y 
other groups dGrise et al.ll2008ll2009Hlmpiombato et alj|201 lh . 



M82 X-l is located near o r within a super star-c luster MCG- 
11 with an age of 7-12 Myr dMcCradv et al.ll2003l) . NGC 1073 
IXO 5 appear s near a ring of star formation of 8-16 Myr old 
dKaaretll2005l) . N GC 1313 X-2 is found near a star cluster with 
an age of 20 Myr dLiu et alJ2007tlGrise et alJl2008l) . IC 342 X- 
1 lies in a fiel d with relatively old s tellar objects, mostly older 
than 10 Myr dFeng & Kaaretil2008l) . The stellar field around 
NGC 1 3 1 3 X- 1 is at least 30 Myr old dYanget alJ201ll) . Holm- 
berg IX X-l is suggested to be pa rt of a loose cluster with an 



age < 20 Myr dGrise et al.ll201ll) . All those findings suggest 



that ULXs are in relatively old stellar environments with ages 
of at least 10 Myr. 

An exceptional c ase is provided by a transient ULX recently 
discovered in M83 ( Soria et al ]|2010ah . It is located in an inter- 
arm region with a stellar population older than ~ 10 8 yr. HST 
observations when the system was probably in quiescence con- 
strain the donor star to be < 5M . Instead, Gemini observations 
during the recent outburst revealed a bright, blue counterpart 
with absolute B magnitude ~ -4, due to the irradiated outer 
disk and/or irradiated donor star. 

Identification of the optical counterparts to UL Xs is the first 



step fo r a potential dynamical mass measurement. Pak ull et al 



d2006l) found that the broad He 11 ,14686 line from NGC 1313 
X-2 changed its velocity by about 380 km s" 1 between two ob- 
servations. If this is due to binary motion of the disk, an upper 
limit on the BH mass of 50 M can be placed. There have 
been a few attempts to obtain the mass function of a few ULXs 
using emission-line velocity curves. However, the lines were 
found to vary randomly instead of showing an ellipsoidal mod- 
ulation; thu s, they did not provi de dynamical constraints on the 
BH masses (Roberts et al HoT It). The orbital period is a key pa- 
rameter for dynamical mass measurement. With current facili- 
ties, it may be more feasible to search for orbital periods in X- 
rays and then confirm them in optical for extragalactic sources, 
like the cases of M33 X-7 and IC 10 X-l. 
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7. Jets, outflows and bubbles 

7.1. X-ray photoionized nebulae 

We naturally expect ULXs to ionize the surrounding in- 
terste llar medium (see Tarter et aDll969i: Kallman & McCrav 
1982) for the basic principles of X-ray ionization. The main 



difference between UV ionized H n regions and X-ray ion- 
ized nebulae is the lack, in the latter, of a sharp boundary 
(Stromgren sphere) bet ween the ionized and neutral plasma 
(Pak ull & Mirionil 120021) . Instead, there is an extended warm 
region where weakly ionized atoms coexist with collisionally 
excited neutral species. As a result, an X-ray ionized nebula 
is characterized by high-ionization emission lines such as He n 
/14686 and [O iv] at 25.89 fim in the inner regions, and low- 
ionization forbidden lines such as [O i] /16300 in the outer re- 
gions. In particular, He n /14686 acts as a p hoton counter of the 
ionizi ng flux between 54 eV and « 200 eV jPakull & A ngebault 
1986). By comparing the true source flux inferred from the 



He ii ,14686 line with the apparent flux measured at energies 
> 300 eV with Chandra or XMM-Newton, we can test whether 
the X-ray emission is isotropic or beamed. 

The N159F nebula, powered by the sMBH LMC X-l, is 
the best-known example of a large -scale (~ 10 pc) X-ray ion - 



ized nebula in Galactic system s ( Pakull & Angebault 1986b . 
IPakull & Mirionil d2002i 120031) opened a new field with their 
systematic search for similar nebulae around ULXs. They 
found a spectacular case around the ULX in M8 1 group dwarf 
galaxy Holmberg II. The Holmberg II X-l ionized nebula (now 
known as "Foot N ebula") has a He n /14686 luminosity « 2.5 - 
2.7 x 10 36 erg s" 1 dPakull & Mirionill2002t Ikaaret et alj|2004h . 
30 times higher than the LMC X-l nebula, which implies an 
average^ ionizing X-ray flux ~ 4-6 x 10 39 erg s _1 . This re- 
sult has been confirmed by further optical/infrared spectro 



scopic stu dies (Lehma nn et al.l l2005t lAbolmasov et al.l 12007 



Berg hea et al.1 120101) . On the other hand, the apparent X-ray 
luminosity of Holmberg II X-l has been seen to vary over 
the years between ~ a few 10 39 erg s" 1 a n d « 3 x 10 40 erg 
s" 1 dZezas etaD 119991: iMivaii et all l200ll iGoad etaD 12006 

aiJ boia r 



Grise et ; 



Caballero-Garcfa & Fabian 2010). The con 



sistency between the two independent measurements gives the 
strongest evidence to date that beaming is negligible (l/b < 
3) for this ULX; this favors models with relatively large BH 
masses (Section 0. Similar arguments have also been applied 
to the photoionization-dominated nebula around NGC 5408 X- 
1 . From the flux in high-ionization fines su ch as He n /14686 
and [Ne v] ,13426, iKaaret & Corbell d2009l) inferred a lower 



limit on the isotropic X-ray luminosity *3x 10 erg s , only 
a factor of 3 lower than the directly measured X-ray luminos- 
ity. As in Holmberg II X-l, this suggests that beaming factor 
is < 3. Mo re examples of X-ray ionized nebulae in ULXs ar e 
repo rted bvlPakull & Mirionil d2002l) . lAbolmasov et all d2007l) . 
and lKaaret et al.1 d2010l) . although in those cases the interpreta- 
tion is more model dependent, because the observed line spectra 
are a mix of photo-ionized and shock-ionized plasma emission. 



5 over the recom bination time of He ++ in the nebula, estimated as ~ 3000 yr 
iKaaret et all2004l) . 



7.2. Shock-ionized nebulae 

It is often hard to distinguish between X-ray photoioniza- 
tion and shock excitation, for ULX nebulae — just as it is gen- 
erally difficult to separate the contribution of a starburst nu- 
cleus and/or a low-luminosity AGN in a LINER. Diagnostic 
diagrams with several line ratios are comput ed with shock and 
photoionizatio n codes such as mappings m dAllen et al J 12008) 
and cloudy dFerland et al .119981) . Shock-ionized plasma is char- 
acterized by high ratios of low-excitation lines such as [O 1] 
,16300 [S n] AA6717,6731 and [N 11] ,1,16548,6583 over Hor. 
For shock velocities v s > 80-100 km s _1 , the high-ionization, 
forbidden transitions [O m] ,114959,5007 are strong ([O in] A 
5007/Ha * 1 for v s > 100 km s" 1 ), whil e He 11 ,14686 indi- 
cates v s > 300 km s" 1 ( Dopita et alj 1984]) . The luminosity in 
the Hy6 line is o nly a function of the total s hock power P w and 
shock velocity dDopita & Sutherland! 19951) : for typical speeds 
~ 100-300 km s _1 , P w ~ 100-300L^. 

Supershells were first found around N GC 6946 M F 16 
dBlair&Fesenlll994 and Holmberg IX X-l dMillerll 19951) . and 
were initially speculated to be powered by unusually energetic 
SN remnants. The X-ray emis sion could not be from th e neb- 
ulae, because of its variability (Roberts & Colbert 2003). Fol- 
lowing the pioneering work of Pakul l & Mirionil d2002l bdol . 
several more huge shock-ionized nebulae have been identi- 
fied around ULXs ("ULX bubbles", see Figure [5] for exam- 
ples), with characteristic diameter ~ 200^1-00 pc, expansion 
velocity ~ 100-200 km s , characteristic age ~ 0.5-1 Myr, 



mechanical power ~ 10 —10 erg s (independently deter- 
mined from the bubble expansion speed and from the HB lumi- 

1 11 — 1 — 11 l 

nositv; lWeaver et alJll977| ;|Dopita & Sutherland 1995) an d en- 



ergy content ~ 10 53 e rg dRoberts et alJl2003c Iprise et alj|2 006: 
Ram sey et all l2006t lAbolmasov et alJ 1200 7: Pak ull & Grise 



2008i lKaaret et al.l2010UMoon et alJ201 lllRussell et al.l201 lh . 
The two best-studied examples are t he bubbles around NGC 
1313 X-2 and Holm berg IX X-l dPakull & Mirionil 12002 : 
Pakull & Grisel 120081) . ULX bubbles resemble ordinary SN 



remnants, but are an order of magnitude larger and 10 to 100 
times more luminous. They cannot come from (hypothetical) 
single hypernova explosions: first, because the surrounding 
stellar populations are already too old (typically > 10 7 yr) and 
without O stars; secondly, because it is unlikely that a lower- 
mass binary companion that survived the hypothetical hyper- 
nova explosion could evolve t o fill its Roche lobe an d feed the 
ULX already within < 1 Myr dPakull & Grisell2008l) . Invoking 
10 to 100 ordinary SN explosions over the last 1 Myr would re- 
quire a massive st ar cluster with M > a few 1 5 M o {e.g., from 
STARBURST99; iLeitherer et ai]|2010l 1 19991) . However, none 
of those ULXs is associated with a massive star cluster: at most, 
they are inside ordinary OB associations with age ~ 10-20 Myr 
and total stellar mass ~ 1O 4 M . The only plausible conclusion 
is that ULX bubbles are inflated by a continuous jet or outflow 
from the ULX, with a power comparab le to the radiative output 
dPakull et alJl20()3lPakull & Grise1l2008l) . 

Jet-inflated bubbles around accreting BHs are found in our 
Galaxy, too. The best-known example is the W50 ne bula 



around SS 433 dDubner et all 119981: iGoodall et alj 1201 lh : a 



spherical shell with two lateral extensions, the ears, that are 
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Figure 5: Large Ha emission nebulae aroud a selection of ULXs and in the jet-dominated source NGC7793 S26. Clockwise, from top left: HolmberglX X-l 
(Subaru/FOCAS image from 2003, credit F. Grise); NGC 1313 X-2 (VL77FORS1 image from 2003, credit F. Grise); IC342 X-l (HST/ACS image from 2005); 
NGC 7793 S26 (HS77WFC3 image from 201 1). In each panel, the X-ray source position is marked with a circle of 0.5" radius. 
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inflated by interaction of the jets with the interstellar medium. 
The characteristic diameter along the major axis is a; 100 pc. 
Both the jet power and the radiative luminosity are > 10 39 
erg s _1 . Indeed, it is suggested that SS433 would ap pear as a 



beam ed, super- Eddington ULX if it was seen face-on (Fabrika 
12004 . Another loca l example is the 5-pc bubble inflated b y the 
sMBH CygnusX-1 dGallo et alJIioollRussell et al.ll2007h : the 
shock velocity determined from optical emission line ratios is 
~ 100-360 km s _1 ; the required jet power » 10 37 erg s _1 is 
~ 0.3-1 times the X-ray luminosity. 



Eddington power. Mechanically-dominated, radiatively ineffi- 
cient quasars in the early universe may grow much faster than 
Eddington-limited sources, and would provide strong mechan- 
ical rather than radiative feedback onto the surrounding gas. 
Thus, understanding the inflow/outflow/radiation processes in 
ULXs has wider astrophysical relevance. Moreover, radia- 
tive and mechanical feedback from HMXBs and ULXs them- 
selves may have been an additional, important source of heating 
and reionization of the intergalactic medium at high redshifts 
dMirabel et alJl201lb . 



7.3. Jet-dominated sources 

The ULX bubbles discovered so far have a quasi-spherical 
appearance, with no direct signature of relativistic, collimated 
jets (as in the Galactic source SS 433). The first extra-Galactic 
super-bubble with collimated jets was recently found in th e out- 
skirts of the spiral galaxy NGC 7793 JPakull et al.ll2010b . This 
source (Figure 0, known as S26, has: X-ray and optical hot 
spots where the jet and counter-jet interact with the interstellar 
medium; an optical nebula with a size of « 300 x 150 pc, ex- 
panding at as 250 km s _1 ; an associated X-ray and radio nebula 
with bright radio l obes. The long-term-average jet p ower ~ a 
few 10 40 erg s" 1 dPakull et alJl2010t ISoria etal]l2010bl) . over a 
characteristic age as 2 x 10 5 yrs. The X-ray luminosity of the 
central BH is only 7 x 10 36 erg s _1 . This might be because the 
BH emission is highly obscured and we only see a scattered 
component; or because it is beamed away from us; or because 
the BH is a transient, currently in a radiative-faint state, but may 
appear as a ULX at other times. Whatever the reason for the 
current radiative dimness, it is important to note that NGC 7793 
S26 has had collimated jets with a long-term average power 
at least of order of the Eddington power (and probably much 
greater) of a stellar BH (sMBH or MsBH). In other words, it 
is a jet-dominated system at super-Eddington mass accretion 
rates. There is no physical principle that prevents mechanical 
power to exceed the Eddington luminosity. However, the fact 
that it is collima ted in a large-scale jet may be surprising. The 
canonical view dFender et al. 2004) of BH accretion states — 
based mostly on the behavior of Galactic transients — is that col- 
limated, steady jets are found at m < 0.03, when accretion is ra- 
diatively inefficient (low/hard state). Whether and how steady, 
collimated jets can be launched during super-Eddington accre- 
tion phases remains an unsolved theoretical problem. A pos- 
sible solution c omes from the radi ation-magnetohydrodynamic 
simulations of Take uchi et al.l(12010l) . who find that a jet can be 
accelerated up to mildly relativistic speed by radiation pressure, 
and collimated by the Lorentz force of a magnetic tower made 
of wound-up toroidal magnetic field lines. The magnetic struc- 
ture is itself kept collimated by the geometrically thick accre- 
tion flow that dominates in the innermost region when m > 1 . 

An ultrapowerful (rather than ultraluminous) system such 
as NGC 7793 S26 is the local-universe analog of a recently 
discovered class of q uasars, dominated by mechanical power 
(Punslv 2007, 1201 II) . In those quasars, the jet power is up 
to 25 times higher than the bolometric radiative luminosity; 
in some of those systems, the jet power is near or above the 



7.4. Radio bubbles 



Some ULX bubbles (from both the photoionized and shock- 
ionized variety) are associated to extended radi o emission. 
The b est-studied example s are Holmberg II X-l dMiller et al ' 
200 5| and NGC 5 408 X-l dKaaret et alj2003l:ISoria et alj20Qi 



Lang et al. I l2007h . Both radio nebulae are dominated by op- 



tically thin synchrotron emission, somewhat similar but much 
larger (« 60 x 40 pc and « 40 pc, respectively) and more lumi- 
nous than typical radio SN remnants. They are also an order of 
magnitude more lum inous than the radio nebula around SS 433 
(Dub ner et al.l ll998). which well matches the relative values of 
radiative and mechanical core power in those sources. The total 
power carried by synchrotron-emitting electrons has been esti- 
mated to be between ~ a few 10 49 erg s _1 and « 10 51 erg s _1 , 
depending on the very uncertain lower-energy cut-off of the dis- 
tribution. In any case, both nebulae are much more energetic 
than W50 around SS433 (~ a few 10 46 erg s -1 : bubner et all 



1998) or than the extended lobes of typical Galact ic micro- 
quasa rs such as GRS 1758-258 (» 2 x 10 45 erg s"': |Hardcastle 
| 2005l) NGC 7793 S26 also has a large radio nebula dSoria et al. 
BOlObh . interpreted as a combination of flat-spectrum free-free 
emission (mostly in the central region) and steep-spectrum syn- 
chrotron emission from the radio hot spots and lobes. The en- 
ergy stored in relativistic electrons is ~ a few 10 50 erg s _1 . 

The energy content in relativistic electrons for typical ULX 
nebulae is only ~ 0.01-0.001 of the core power — in particular, 
of the mechanical power when it can be reliably estimated from 
bubble size or optical line fluxes — integrated over the charac- 
teristic age of the bubble. The main reason is that most of the 
input energy is transferred to non-relativistic elec trons, protons 
and nuclei dPakull et alJbOlOtlSoria et alj|2010bl) . The fraction 
of energy in non-radiating particles to that in relativistic elec- 
trons in a system such as NGC 77 93 S26 is similar to the values 
found by (Cavagn olo et alJEoiol) for a sample of AGN, Seyfert 
and giant elliptical galaxies. For the same jet power, AGN, 
ULXs and Galactic microquasars may be radio louder or radio 
quieter depending on the amount of plasma entrained by the 
jets. We speculate that an interesting test may come from com- 
parisons of the radio loudness and lobe energy in low/hard state 
micr oquasars, probably dom inated by magnetically accelerated 
jets dKomissarov et al. 2007), and in ULXs, poss ibly dominated 
by ra diation-pressure accelerated jets/outflows dTakeuchi et al 
2010), with more entrainment. 
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8. Summary and future prospects 

There is no doubt that our understanding of the phenomenol- 
ogy of ULXs has greatly advanced in the past decade thanks 
to Chandra, XMM-Newton, and many other space- and ground- 
based facilities. However, the key physical issue is still unre- 
solved: that is, whether they are powered by IMBHs or nor- 
mal stellar BHs. A few sources have emerged as strong IMBH 
candidates: in particular, ES0243-49 HLX-1 and M82 X- 
1 . On the other hand, for the majority of ULXs, there is no 
strong theoretical need nor compelling observational evidence 
for IMBHs. In the absence of direct mass-function measure- 
ments from phase-resolved optical spectroscopy, we still have 
to rely on X-ray spectral and timing modeling and other indi- 
rect clues. Reviewing the results that we have collected here, 
we suggest the following conservative statement: 

ULXs are a diverse population; MsBHs with moder- 
ate super-Eddington accretion seem to be the easiest 
solution to account for most sources up to luminosi- 
ties ~ a few 10 40 erg s~ l ; strong beaming (\/b> 10) 
can be ruled out for the majority of ULXs; IMBHs are 
preferred in a few exceptional cases. 

Some of the following kinds of observations may provide fur- 
ther breakthroughs: 

a) measurements of the fast timing behavior down to a time 
scale of sub-seconds may probe the inner part of the inflow: 
ideally, we would search for possible high frequency fea- 
tures (breaks and QPOs) that are found in Galactic BHs at 
frequencies ~ 10 2 Hz. This would require a new generation 
of X-ray telescopes with large effective area. 

b) determining the relative contribution of thermal emission 
and Comptonization component is a key test for competing 
accretion models. For that, multilayer-coated X-ray tele- 
scopes with good sensitivity up to a few tens of keV are 
needed. 

c) detection of, or a more solid upper limit to broad iron K lines 
may constrain the properties of the inner disk. 

d) High-resolution X-ray spectroscopy with gratings or mi- 
crocalorimeters may constrain the metal abundance and ion- 
ization state of the inflow, the mass and energy of the out- 
flow, and probe the absorption from the surrounding inter- 
stellar medium. 

e) the ULX sample within ~ 15 Mpc is almost complete, and 
it does not contain enough sources to permit statistical stud- 
ies of the luminosity distribution above « a few xlO 40 erg 
s _1 . A homogeneous all-sky survey, complete for ULXs up 
to ~ 50 Mpc would allow us to determine whether or not 
the luminosity distribution extends unbroken up to at least 
10 41 erg s _1 . Such survey will hopefully be successfully 
conducted by eROSITA; 

f) repeated snapshot monitoring of nearby ULXs over the next 
few years will extend their long-term lightcurves enough to 
prove that their duty cycles and transient behavior are dif- 
ferent from those of Galactic BHs. This can be easily done 
with smaller X-ray missions; 



g) phase-resolved spectroscopy would lead to the measurement 
of radial velocity curves and mass functions. Attempts to do 
so with 8-m class telescopes have so far been inconclusive 
but have suggested that it will be easily done with future 
30-m class telescopes; 

h) integral-field spectroscopic studies of star forming galaxies, 
as well as low-frequency radio studies, may reveal fainter, 
older ULX bubbles or lobes, still expanding into the inter- 
stellar medium but no longer energized by a central source. 
This would constrain the duration of the ULX active phase; 

i) searching for compact radio jets down to a sensitivity of 
~ 1//Jy will constrain the BH mass, in relation to other 
classes of accreting BHs in the fundamental plane. This 
may be done in the near future with the eVLA and with SKA 
Pathfinders, and in the more distant future with the SKA it- 
self. 

In conclusion, we gather from this list that there is more than 
one way to skin a cat. Many fundamental aspects of ULX 
physics will be understood with instruments already planned 
or under construction. Thus, we anticipate a fruitful decade of 
ULX studies beyond the Chandra and XMM-Newton era. 
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